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Abstract for the project 


Computer simulation of the condensation and sublimation of 
CO2 frost in the martian polar caps has been fairly successful 
in reproducing the annual cycle in atmospheric pressure observed 
by the Viking Landers. However, these studies have not been 
able to uniquely explain why CO2 frost survives southern 
summer. We propose to study several mechanisms which might 
explain this phenomenon: an improved treatment of cloud and 

atmospheric-dust effects; the decrease in the importance of heat 
conduction due to existing CO2 ice; an improved treatment of 
sublimation microphysics; the variance in ice albedo with solar 
zenith angle, ice age, and dust/ice ratio; the bidirectional 
reflectance of ice; snowfall accumulation in addition to surface 
frost formation; insulating residues on the ice; wind shifting 
of ice and brightening of ice under exposure to sunlight. To 
perform this study, we will integrate two existing models of the 
martian polar caps and combine them with an accurate model of 
martian atmospheric radiation and parameter izat ions of the 
processes described above. The model will be used to identify 
combinations of parameters which will allow a reproduction of 
seasonal observations of atmospheric pressure and the 
latitudinal and seasonal extent of the polar cap, and the 
hemispherical asymmetry in the residual polar caps. A better 
understanding of the martian polar caps will have direct 
implications for our understanding of the climate and atmospher- 
ic dynamics of Mars and the Earth, as well as provide insight 
for the observations planned for the Mars Observer mission. 
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II . 


Summary of Research to date 


II. 1. Jakoskv Model Acquired 

The Principal Investigator, Dr. Lindner, flew to Colorado 
in September 1989 to meet with the Co-Investigator, Dr. Jako- 
sky. The polar cap model of Dr. Jakosky was obtained. Dr. 
Lindner became familiar with the code in discussions with Dr. 
Jakosky, and through exercising the model on Dr. Jakosky' s 
computer system. Dr. Lindner and Dr. Jakosky also held in-depth 
discussions of the overall approach to be taken on the project 
to meet the objectives. 

II. 2. Lindner and Jakoskv Models Integrate d. 

The Jakosky model was loaded onto the AER computer system 
and combined with the Lindner model, as discussed in the 
proposal. The hybrid model was debugged of errors. 

II. 3. Comparison to Kieffer's su rface model 

Kieffer (1990) described a model which is similar to ours, 
which he used for ice metamorphism studies. I set up the code 
to use the same inputs as Kieffer did in his Figure 3, and 
plotted the results on the same figure (deleting all but the 
case I compared to). The agreement is good, as can be seen in 
Figure 1. Since I don't know exactly what all of his input 
parameters were (such as emissivity, seasonal ice albedo, etc.), 
I consider the agreement a validation of our code. The only 
different behavior is why his figure doesn't show a seasonal 
temperature "pulse" which makes it's way down in depth (i.e. 
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1. Temperatures within the north polar ice cap (shown are 
the annual minimum and maximum). Lindner and Jakosky 
values are for the current study (L — 260 and 280 
degrees), and Kieffer values are taken from Kieffer 
(1990). The albedo of the residual ice is 0.41, and the 
thermal inertia is 0.04 cal cm-2 sec-0.5 K-l. Surface 
values are shown as 1 cm depth. 



when the surface is at a maximum, the temperature at 10 m should 


be at a relative minimum, as our results show). Perhaps Kieffer 
is showing the maximum temperature for all year at each depth, 
rather than a snapshot in time when the surface reaches a 
maximum (as I've plotted). 

II. 4. Ice Microphvs ics Parameterization 

Section 2.3 of the proposal discussed the inclusion of 
improved ice microphysics in the model. Discussions were held 
with Dr. Gary Clow as to the nature of his ice microphysics 
models and the results he has obtained recently. In the 
proposal, we mentioned that the depth to which radiation pene- 
trates Into the cap has never been included in polar cap 
sublimation models, although it could be Important. We have 
done calculations that show that this is not important for thick 
C0 2 ice, or thin C0 2 ice overlying a dust surface. For the 
thick ice case, the Integral of sublimation Is Independent of 
where it occurs within the ice. The distribution of absorption 
within the polar cap is important for studies of the Internal 
physics of the cap, but not for overall sublimation rate. For 
the thin C0 2 ice over dust case, the surface thermally 
reradiates the absorbed solar radiation primarily back to the 
polar cap where most is reabsorbed. The top few grains of dirt 
absorb the incoming solar radiation which penetrates the polar 
cap, and it is much easier for these top few grains to conduct 
or radiate heat back to the overlying polar cap than deep into 
the surface . 
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A thin CC >2 ice polar cap overlying a residual water ice 
deposit is different. Solar radiation which passes through the 
CO 2 ice cap may penetrate quite deeply into the 1^0 ice. This 
will effect C0 2 ice lifetimes over the residual polar cap, but 
won't affect overall polar-cap recession or the annual variation 
in surface pressure, both of which current models reproduce. I 
incorporated this effect in the model by apportioning the 
absorbed solar flux between the sublimation of the ice, and the 
heating of the subsurface. Clow (1987) has already calculated 
the absorption of solar flux as a function of depth in C02 ice 
on Mars, and I have reproduced a figure from his paper in Figure 
2, which we used in the model. Based on Clow's figure, I 
obtained a % of solar flux absorption as a function of depth for 
41 depth steps. The code already calculates the thickness of 
C0 2 ice, and I use this depth with Clow's work to determine how 
much sunlight penetrates the seasonal ice. This fraction that 
penetrates the seasonal ice is further apportioned between the 
individual subsurface layers, again using Clow's figure. 

Both Kieffer (1990) and Moore (1988) have deduced that ice 
grain radii in the residual cap are over 100 pm and that dust 
concentrations are less than 1/1000, so any one of Clow's 
profiles could be legitimate. I've tried them all, and found 
the largest deviation was for clean, fine ice. This is 
understandable, since solar radiation penetrates the deepest in 
clean, fine ice. In our work, we have assumed that the seasonal 
ice was fine, and the residual ice was coarse. 
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Fig. 3. The integral of net shortwave flux with re- 
spect to depth for clean and dusty snow, using the 
mean-annual incident solar flux at latitude —38°. 
Mean-annual temperature profiles in snow are propor- 
tional to the flux integral / s (z), which is generally 
greater for coarse-grained than for fine-grained snows. 
Curves (A) correspond to a snow with ice grain radii of 
50 /xm and bulk density 50 kg m -3 while curves (B) are 
for a snow with 1000-am ice grains and a density of 400 
kg m -3 . {jProm Clow (1987^ 
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I have plotted in Figure 3 the difference in the depth 
profile of temperature between the code with no penetration of 
light and the modified code using penetration appropriate for 
clean ice for Lg - 260° (just before the last of the seasonal 
ice sublimes). The radiation which penetrates the cap heats the 
subsurface by up to 3°K. Most of the radiation gets absorbed 
near the surface, but that easily conducts away to the 
surface. The greatest heating occurs at 1 m depth. The scale 
on the right side of the figure is for the "difference” curve. 
The surface temperature is plotted as 1 cm depth on this plot. 

Figure 4 shows the profiles for Lg«280° , when the maximum 
surface temperature occurs, as well as the maximum deviation at 
depth between the codes (9K). The surface is actually a couple 
of degrees cooler for the penetrating radiation case, as we 
expected. Again, the maximum subsurface heating occurs at 1 m 
depth . 

Figure 5 shows how the surface temperature varies with Lg 
for both codes. Surface temperatures are cooler from Lg - 270° 
to Lg - 310° for the penetration code, since solar radiation is 
absorbed at depth, and must conduct to the surface. From 
Lg - 310° to 360° , the surface is actually warmer than for the 
original code, due to the increased conduction by the warmer 
subsurface . 

Figure 6 shows depth profiles of temperature for every 10 
sols through the year for clean ice penetration, and Figure 7 
shows the same for the no-penetration code. Again, the heating 
of the subsurface can be seen for summer (just before the 
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Year 6; Solar Longitude: 260 Degrees 
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3 ' T® m P erature versus depth in the south polar cap for L 
of 260 degrees simulated including the effect of light 
penetration for "clean" ice, and simulated without 
penetration. The difference between with and without is 
also plotted, using the scale on the right side. Surface 
va -*“ ues sre shown at 1 cm depth. 
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Fig. 4. Temperature versus depth in the south polar cap for L 
of 280 degrees simulated including the effect of light 
penetration for "clean" ice, and simulated without 
penetration. The difference between with and without is 
also plotted, using the scale on the right side. Surface 
values are shown at 1 cm depth. 
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Fig. 5. Surface temperature of the south pole versus Solar 
Longitude simulated including the effect of light 
penetration for "clean" ice, and simulated without 
penetration. The difference between with and without is 
also plotted, using the scale on the right side. 
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Fig. 6. Temperature versus depth in the south polar cap for 

f every 10 sols through one year simulated including the 

effect of light penetration for "clean" ice. Surface 
values are shown at 1 cm depth. 
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Fig. 7. Temperature versus depth in the south polar cap 
every 10 sols through one year simulated without the 
light penetration. Surface values are shown at 1 
depth. 
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surface ice sublimes away) and in the warm subsurface at the 
— maximum temperatures. We also see that very little annual 

variation in temperature occurs at 4 m depth. 

Figures 8 and 9 show the annual variation in temperature, 
with each curve representing a different depth (at z-0, 2.8, 

9.5, 19.0, 32.2, 50.6, 76.5, 112.7, 163.4, 234.4, 333.7, 472.8, 

_ 667.5, 940.1, 1321.8, 1856.2, 2604.2, and 3651.5 cm). The 

deepest depths show the least variation (almost a straight line 
' at 36 m. The ice-covered heating of the subsurface can be seen 
between Lg - 180° and 270*. 

Figure 10 shows the effect of light penetration on surface 
mm frost between the two models in the first year of iteration 

(starting from the same initial conditions of an isothermal 
^ surface of 154K, and no frost). Since the sun doesn't shine 

from L s -0° to 180° , no difference is seen then. As the cap 
thins in the spring, the effect of sunlight penetrating the cap 
^ becomes more important (as expected from Clow's figure). Thus 

the ice in the model which includes penetration lasts longer in 
Year 1. 

Figure 11 shows the same for Year 6, by which time the 
models have converged. Now the model which includes penetration 
has less ice all year. This is because the subsurface heats up 
more in the summer due to light penetration, and this heat 
decreases the amount of condensing frost in the fall and winter, 
resulting in an increasing difference between the models from 
L s -0 to 180°. After Lg-180 0 , this difference decreases because 
of light penetration as we saw for Year 1, but not enough to 
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Fig. 8. Surface temperature of the south pole versus Solar 
Longitude, simulated including the effect of light 
penetration for "clean" ice. Each curve represents another 
layer in the model, ranging from the surface to 3.5m depth. 
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Fig. 9. Surface temperature of the south pole versus Solar 
Longitude, simulated without the effect of light 
penetration for "clean" ice. Each curve represents another 
layer in the model, ranging from the surface to 3.5m depth. 
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10. Frost amount at the south pole versus solar longitude 
simulated including the effect of light penetration for 
"clean" ice, and simulated without penetration. The 
difference between with and without is also plotted, using 
the scale on the right side. Results are for the first 
year of model iteration, starting from an isothermal 154K 
surface and lOg cm-3 of frost. 
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Fig. 11. Frost amount at the south pole versus solar longitude 
simulated including the effect of light penetration for 
"clean" ice, and simulated without penetration. The 
difference between with and without is also plotted, using 
the scale on the right side. Results are for the sixth 
year of model iteration. 
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compensate completely, and the Ice sublimes away earlier for the 
model which includes penetration. 

Figures 12 and 13 show the depth profiles of temperature 
for the dirtv Ice case versus the no - pene t r a t ion model. The 
differences in temperature exhibit the same behavior as with 
clean ice, but of a much lower magnitude (as we expected). 

Figure 14 shows surface temperature as a function of Lg for 
dirty ice and the original code, again exhibiting the same 
behavior as for clean ice, but with a smaller magnitude. 

Figure 15 shows the effect of light penetration on the frost 
budget for dirty ice, again showing the same behavior as for 
clean ice, but with a smaller magnitude. 

I have checked to make sure that these deviations in 
temperature are not caused by the modifications to the code. To 
check that these changes weren't altering the code in its 
original form, I ran one case where all solar radiation is 
absorbed in the top mm of Ice, rather than distributed over the 
top meter or so. This case gave the same results as Jakosky's 
original code, and Kief fer ' s model. 

Obviously, the code had to be run with very thin layers at 
the surface (a 5 cm surface layer thickness for this case). 
Thicker layering caused an error, which was determined by 
comparing to the original Jakosky code. Even a 26 cm top layer 
caused a 3 sol increase in the number of ice-free days, and a 
100 cm top layer caused a 15 sol error. This is obvious, since 
the layering needs to be on the same scale as the penetration 
depth. 
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Fig. 13. Temperature versus depth in the south polar cap for L 
of 280 degrees simulated including the effect of light 
penetration for "dirty" ice, and simulated without 
penetration. The difference between with and without is 
also plotted, using the scale on the right side. Surface 
values are shown at 1 cm depth. 
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14. Surface temperature of the south pole versus solar 
longitude, simulated including the effect of light 
penetration for "dirty" ice, and simulated without 
penetration. The difference between with and without is 
also plotted, using the scale on the right side. 
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The results do depend on the input parameters (albedo of 
ice, inertia, albedo of surface, emissivity, etc.). For 
example, a low thermal inertia (0.01) yields a much hotter 
subsurface and dramatic reduction in ice-covered sols due to 
light penetration. Also, the albedo needed to keep ice year- 
round is -0.04 higher for frost emissivity of 0.95 versus 1.0 
(obviously, since a reduction in emissivity reduces amount of 
heat lost, and a higher albedo reduces the absorbed solar 
radiation to maintain balance). The results do not depend 
strongly on which hemisphere is studied. However, the results 
do depend on the type of ice (coarse/fine and clean/dirty). A 
dirty, coarse ice has much less light penetration than a clean, 
fine ice. But all cases result in fairly small changes, and can 
be neglected given the uncertainties in other parameters. 

I did run one case with soil as the underlying surface. As 
we had stated earlier, there is no difference between the no- 
penetration code and a code which allows light to penetrate the 
seasonal cap and get absorbed in the top mm of soil. Almost all 
the penetrating radiation returns to sublime the seasonal ice 
from this depth. 

However, this does imply that the seasonal CO 2 frost 
sublimes away more quickly over the residual cap in the spring 
than over the neighboring soil, and condenses more slowly over 
the residual cap in the fall than over the surrounding soil. 

This might tend to move the residual cap around a bit on a 
climatic timescale, but would probably be a small effect. 
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II , 5, 


Solar zenith angle dependence of albedo, . 


The albedo of ice depends on solar zenith angle (SZA). 

This SZA dependence has been theoretically predicted for C0 2 
frost on Mars by Warren et al . (1990) . I assumed that the co- 

albedo (1-albedo) varied with the SZA dependence given by 
Warren, no matter what the albedo was. This will give 
reasonable albedos when near the 0,7 base albedo used in 
Warren's calculations, but is speculative for albedos far from 
this value. Fortunately, the values considered here are very 
close to this value. 

As mentioned by Warren et al. (1990), his calculations are 
for a flat surface. Because of surface roughness on a natural 
snow surface, the effective zenith angle is usually much less, 
and rarely over 80° (Warren et al . , 1990). I have simply 

adjusted the SZA to be 10° less than given by orbital mechanics 
as a simple approximation for surface roughness. 

I find that the influence of the sun on ice stability is 
significantly reduced in spring and fall, when SZA 's are 
large. This significantly lowers the albedo required to 
maintain ice on the south pole year-round, by about 5%. 

II. 6, All parameterizations considered together, w 1 th the — best 

fit to the data. 

The results discussed up to now dealt strictly with Ice at 
90 ° N and 90°S latitude. The problem with polar cap models to 
date has been that they cannot maintain C0 2 ice on the south 
pole without also maintaining it on the north pole, unless a 


27 


hemispherical asymmetry in ice albedo is assumed. Having noted 
that the effects of clouds and dust can correct this problem 
(Lindner, 1990), and having seen that the SZA effect on albedo 
can allow for reasonable albedos to accomplish this, 1 set out 
to run the model at all latitudes and seasons to see if the 
observed data on polar cap regression and atmospheric pressure 
could be reproduced. 

The Lindner (1985) polar cap model was modified to allow 
for light penetration into the cap, to allow ice albedo to vary 
with age, latitude, hemisphere, dust content and solar zenith 
angle, to include the radiative effects of clouds and dust, to 
allow for diurnal variability, and to include heat transport as 
represented by a thermal wind, as discussed in the proposal. 

Ice albedo was chosen to have either the albedo for old ice or 
new Ice, depending upon whether frost was condensing or 
subliming, with a gradual transition from the albedo of new Ice 
to that of old ice over a period of 10 days. The albedo of old 
ice was usually assumed to be less than that of new ice, 
although several cases of brightening of ice with age of the 
type postulated by Paige (1985) were investigated. Ice albedo 
also was allowed to vary linearly with distance from the edge of 
the polar cap, being up to 10% higher if the edge of the polar 
cap were 50° in latitude away. The transport of dust into the 
winter polar region is ineffective and there are no sources of 
dust when the soil is ice-covered. Changes in albedo due to 
dust content within the ice and hemispherical asymmetries in 
albedo were allowed to vary at random to obtain the best fit to 
the data . 
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The radiative effects of clouds and background dust (i.e. 
other than during dust storms) as computed by Lindner (1990) 
were incorporated. The cloud opacity was assumed to be 
hemispher ically asymmetrical, with an opacity of 0.5 in the 
north and 0.2 in the south, and to vary linearly with the annual 
cycle in atmospheric surface pressure observed by the Viking 
Lander. Clouds were assumed to exist only where CO 2 frost 
existed. The annual cycle in background dust opacity was taken 
from Pollack et al. (1977). Dust opacity was assumed to vary 
linearly with distance from the edge of the polar cap, so that 
it was only half that at the equator if the edge of the polar 
cap was 50° in latitude away. Both background dust and cloud 
opacity were also assumed to vary linearly with the changes in 
airmass due to changes in elevation. The radiative effects of 
CO 2 vapor were also included (NIR absorption of solar radiation 
and emission of IR radiation). The scientific basis for all of 
these approximations were discussed in Lindner (1990) and in the 
proposal to this contract and will not be repeated here. 

Heat transport was simply incorporated by the use of the 
thermal wind approximation. Future versions of this model will 
use a full GCM calculation. The thermal wind is based on the 
temperature difference across the edge of the polar cap, and 
various efficiencies for heat transport were chosen for model 
comparisons . 

The model was also run with time steps of 1/50 of a day, 
small enough to account for diurnal variations in solar zenith 
angle, surface temperature, and frost condensat ion/sub lima- 
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tion. The model was run for every 2° in latitude to account for 
latitudinal variation in model parameters. An albedo of 0.45, a 
thermal Inertia of 0.03 cal cm' 2 s' 1 / 2 K' 1 and a surface density 
of 0.93 g ci"^ were chosen for the residual polar caps (Kieffer, 
1990; Jakosky and Haberle, 1990), while values for bare soil of 
0.25 for albedo, 0.006 cal cm* 2 s' 1 / 2 K' 1 for thermal inertia, 
and 1.5 g cm"^ for surface density were chosen for other 
latitudes. Frost temperature was allowed to vary with 
atmospheric pressure. 

I ran hundreds of combinations and permutations of model 
parameters. Rather than present hundreds of graphs, charts, 
etc., allow me to summarize the results and show the best fit to 
the data. The penetration of light into the polar cap was not 
deemed to be important, as discussed in section II. 4. The 
division of albedo between old and new ice was also not very 
important. The fits to the data were not markedly different 
whether old ice was chosen to be higher or lower than new ice. 
The dependence of ice albedo on SZA and the radiative effects of 
clouds and background dust were only important right at the pole 
itself, and not for the overall polar cap regression or 
atmospheric pressure, as postulated by Lindner (1990). The 
radiative effects of global dust storms were not investigated, 
and could be important. Running the model including diurnal 
variability versus in a diurna 1 ly - aver aged mode also made little 
difference to the fit to the data. The inclusion of heat 
transport was significant, and therefore warrants the inclusion 
of more accurate GCM simulations. 
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Figures 16 through 19 show the regression of the polar caps 
as predicted by the model versus observations. Agreement is 
good, particularly in the northern hemisphere. An exact 
comparision is difficult, considering that the edge of the polar 
cap is usually patchy and ill-defined, in large part due to 
terrain. Further, observations of the polar cap are conducted 
during the daylight. Near the cap edge, ice frequently forms at 
night and sublimes away completely during the day. Thus, the 
edge of the polar cap is also diurnally- variable . 

Figures 20 and 21 show the comparison of the atmospheric 
pressure predicted by the model for the elevations of the Viking 
Landers to the actual observations. The agreement is good, 
although the annual cycle in pressure predicted by the model is 
slightly out of phase with the data. This is a problem noted 
with all runs of the model, no matter what sets of input 
parameters were used. The model is still ignoring or poorly 
treating some process which forms ice sooner in the fall, and 
sublimes it away more quickly in the spring. Some suggestions 
for what these may be are presented in the section on plans for 
future work. 

Figures 22 and 23 show model predictions for the amount of 
frost at particular latitudes throughout the year. These values 
can be converted to a depth of frost by using a density of 0.93 
g cm . Frost forms as equatorward as 50°S and 54°N latitude 
with typical maximum depths of a meter. Note that frost forms 
later and sublimes earlier at lower latitudes. Note also that 
once ice begins to sublime, it does so rapidly. 


31 


CAP RADIUS 


| I \ I I I I I I I | IT 


T 


T 


I 


▼ 1931 (Zvactfcl.ftt 

0 1*S< (tV4»Aki «t Al. <1*14)1 
O 1171 (Ja*4« 4 Loom* (1)13)] 

▲ 1973 (I*4«4li «t a 1. (1)14)1 
A 1)73 4 Low* (1)13)] 

V 1)77 (Jams 4 Low* (1)13)] 

* 1)44 £1*444*1 At 41.(1)44)1 
— K 44 A (Pltchb«cK«r 4t 4l. (1)0)1 

♦1)44 (1*444*1 4t 4l.(pr4««Dt 
*e4«ar4»«ntft> 1 



0.0 1, I 

220 230 240 250 260 270 280 290 300 310 

SOLAR LONGITUDE 


g. 16. The regression of the south polar cap, as observed for 
various years (taken from Iwasaki et al . , 1990) and as 

simulated by our model (thin line), as a function of the 
aerocentric longitude of the sun (L s ). The cap radius is 
that whcih would be measured on a polar s te re o graph ic 
projection of the south polar region; the units of the 
radius are fractions of the planetary radius of Mars. 
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Fig. 17. The regression of the south polar cap, as observed in 

1986 (solid circles), 1971 (crosses), and 1977 (plus signs) 
[taken from James et al . , 1990] and as simulated by our 

model (thin line), as a function of the aerocentric 
longitude of the sun (L g ). Two-sigma error bars are 
indicated for the 1986 data; the errors are smaller for 
the denser 1971 data and for the 1977 Viking data. 
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Fig. 18. The regression of the north polar cap, as observed for 
various years [taken from Iwasaki et al. , 1982; symbols 

refer to Dollfus (1973), Baum (1974), Br iggs ( 1 9 7 4 ) , James 
(1979), Capen and Parker (1981), Iwasaki et al . (1982)] and 

as simulated by our model (thin line), as a function of the 


aerocentric longitude of the sun (L ). 
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Fig. 19. The regression of the north polar cap, as simulated by 
our model (thin line) and as observed, as a function of the 
aerocentric longitude of the sun (L g ). Data are taken from 
James et al . (1987); also showing data of Iwasaki et al. 

(1979; 1982). The symbols stand for 1977-1978 and 1980 

Viking data (x's and +'s respectively) and ground-based 
1975-1980 data (circles). 
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PRESSURE (mb) 



TIME (sols) 


g. 20. Daily mean surface pressure at Viking Lander 2, showing 
observations (line with high-frequency oscillations) and 
simulations (smooth line) over the course of one martian 
year starting from the vernal equinox. Observations are 
taken from Hess et al . (1980) , and gaps in the data are due 

to irretrievably lost data. 



PRESSURE (mb) 
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21. Daily mean surface pressure at Viking Lander 1, showing 
observations (line with high-frequency oscillations) and 
simulations (smooth line) over the course of one martian 
year starting from the vernal equinox. Observations are 
taken from Hess et al . (1980), and gaps in the data are due 

to irretrievably lost data. 
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Fig. 22. Mass of CC >2 frost per cra^ of surface at every 2 degrees 
of latitude in the northern hemisphere over one martian 
year starting at vernal equinox. The solid line is 90°N 
1 a t i tude . 



FROST AMOUNT (g cm-2) 


SOUTH POLE, ALBEDO 0.72 



90 180 

SOLAR LONGITUDE (degrees) 


Fig. 23. Mass of C0 2 frost per cm 2 of surface at every 2 degrees 
of latitude in the southern hemisphere over one martian 
year starting at vernal equinox. The solid line is 90 S 
1 a t i tude . 


II. 7. Presentations made at conferences 


Several conferences were attended which were supported in 
whole or in part by this contract. Reprints of papers and 
abstracts published at these meetings are included in the 
Appendix . 

1989 DPS Conference 

Dr. Lindner attended the AAS/DPS conference in Providence 
in October, 1989, and presented a paper entitled "The Martian 
Polar Cap: Radiative effects of ozone, clouds, and airborne 

dust", which discussed the theoretical aspects of this current 
work. The work was well received, and Dr. Lindner benefitted 
from discussions with Dr. Jakosky and others at the conference, 
and from listening to other talks. The abstract was published 
in the Bulletin of the American Astronomical Society, Volume 21, 
p. 979 (1989). NASA contract NASW-4444 paid for expenses. 

LPS Conference 

Dr. Lindner attended the Lunar and Planetary Science 
Conference in Houston in March, 1990 and presented a paper 
entitled "Solar and IR radiation near the martian surface: A 

parameterization for CO 2 transmittance" with T. Ackerman, J. 
Pollack, 0 . B . Toon and G.E. Thomas as secondary authors. This 
paper was published in Lunar and Planetary Science XXI, pp. 696- 
697. NASA contract NASW-4444 paid for some expenses (mostly 
just my labor), with the Lunar and Planetary Institute paying 
the rest of the expenses (airfare, rental car, and some food and 
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hotel) with a grant from the Mars Surface and Atmosphere Through 


Time (MSATT ) program at NASA. 

EGS Conference 

Dr. Lindner attended the European Geophysical Society 
Assembly in Copenhagen, Denmark in April, 1990 and presented a 
paper entitled "CO 2 transmittance in the Mars atmosphere: An 

exponent ial - sum fit for use in mult iple- scatter ing models" with 
J. Pollack and T. Ackerman as secondary authors. This abstract 
was published in Annales Geophysicae, NASA contract NASW-4444 
paid all expenses. 

Atmospheric Radiation Conference 

Dr. Lindner attended the seventh conference on atmospheric 
radiation in San Francisco in July and presented a paper 
entitled "The effects of polar clouds and dust on the radiative 
budget of the martian polar cap". This paper was published in 
the conference proceedings by the American Meteorological 
Society. A. D.O.D. contract paid all expenses except for the 
publication and registration costs related to the Mars work as 
Dr. Lindner was there to present other research as well (Multi- 
spectral Cloud Property Retrieval by B.L. Lindner and R.G. 
Isaacs ) . 

Cloud Physics Conference 

Dr. Lindner attended the cloud physics conference in San 
Francisco in July (concurrent with the atmospheric radiation 
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conference) and presented a paper entitled ”Clouds and Chemistry 
on Mars". This paper was published in the conference 
proceedings by the American Meteorological Society. A. D.O.D. 
contract paid all expenses except for the publication and 
registration costs related to the Mars work. 

1990 DPS Conference 

Dr. Lindner attended the AAS/DPS conference in 
Charlottesville, Virginia in October, 1990, and presented a 
paper entitled "Penetration of Light Into the Martian Polar 
Cap: Implications for the Energy Budget" , which discussed the 

results presented in section II. 4. The abstract was published 
In the Bulletin of the American Astronomical Society, Volume 22, 
p. 1060 (1989). NASA contract NASW-4444 paid for expenses. 

II. 7 Journal Publications 

The text of the proposal for this contract was editted and 
published in the special Mars Polar Processes edition of the 
Journal of Geophysical Research . Solid Earth and Planets, Volume 
95, pages 1367-1379, February 1990. A reprint is included in 
the Appendix . 

A paper by Dr. Lindner entitled " Ozone heating In the 
martian atmosphere " has been submitted to Icarus. A preprint 
will be included in the final report. 
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III. Program of Research for the Next 6 Months 

111.1. Research Tasks 

Seeing that the primary contract objectives have been 
satisfied, and that all contract funds have been spent, the next 
and last 6 months of this contract will be spent preparing the 
results of this contract for publications and presentations, as 
outlined below. 

111. 2. Conferences 

Dr. Lindner will attend the International Union of Geodesy 
and Geophysics Assembly in Vienna, Austria, in August 1991 to 
present a paper entitled ” Mars seasonal Ct^-ice lifetimes and 
the angular dependence of albedo 11 in a special Mars climate 
session. The abstract will appear in the conference 
proceedings . 

Dr. Lindner will also attend the International Symposium on 
the chemistry and physics of ice, held in Sapporo Japan in Sept. 
1991, to present a paper entitled " Why is the north polar cap 
on Mars different than the south polar cap? " in the 
extraterrestrial Ice session. The abstract will appear in the 
conference proceedings. 

111. 3. Publications in Progress 

We are In the process of submitting a short paper to 
Geophysical Research Letters which will present the research 
presented in sections II. 3 and II. 4 of this report. 
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We are also preparing a 
Nature or Science describing 


manuscript to submit to 
the results of sections 


either 
II. 5 and 


II . 6 . 


III. 4 Proposal to continue this work. 

We intend to submit a proposal to 
work. We would integrate more accurate 
transport into our polar-cap model. We 
some processes touched on In the Initia 
which haven't been investigated, such a 
Warren et al. (1990) predictions of ice 
wavelength dependence and bidirectional 
(As predicted by several peer reviewers 
for this contract, I did not request su 
research) . 


NASA to continue this 
simulations of heat 
would also investiga 
1 proposal to this wo 
s terrain, snowfall, 
eraissivity, and the 
nature of ice albedo 
who read the proposa 
fficient funds for th 


t e 
rk 


1 

i s 
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Publications under this contract 


(reprints are in the appendix) 

Lindner, B. L., The martian polar cap: Radiative effects of 

ozone, clouds, and airborne dust (abstract), gull T AjjU. 

As tron . Soc, , 2 1 . 979, 1989. 

Lindner, B. L., Clouds and chemistry on Mars, in Preprints , 
conference on cloud physics , pp . 226-227, held in San 

Francisco, July 23-27, published by Amer. Meteorol. Soc. , 
Boston , 1990a . 

Lindner, B. L., The effect of polar clouds and dust on the 

radiative budget of the martian polar cap, in Preprints r 
seventh conference on atmospheric radiation , pp . 200-201, 

held in San Francisco, July 23-27, published by Amer. 
Meteorol. Soc., Boston, 1990b. 

Lindner, B. L., The martian polar cap: Radiative effects of 

ozone, clouds, and airborne dust, J. Geophvs. R e s, . 9_5, 

1367-1379, 1990c. 

Lindner, B. L., The thermal and frost budgets of the martian 
seasonal polar caps including infiltration of solar 
radiation, submitted to Geophvs. Res. Lett. . 1991a. 

Lindner, B. L., Ozone heating in the martian atmosphere, 
submitted to I c arus . 1991b. 

Lindner, B. L., Why is the north polar cap on Mars different 

than the south polar cap?, Proceedings . International Symp..., 

on the Phvs . and Chem. of Ice , in Press, 1991c. 

Lindner, B. L., Mars seasonal C02-ice lifetimes and the angular 
dependence of albedo, in Proceedings. International Union 

of Geodesy and Geophvs . General Assembly , in Press, 1991d. 
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Lindner, B. L., The polar caps of Mars, in preparation, 1991e. 

Lindner, B. L. and B. M. Jakosky, Penetration of light into the 
martian polar cap: Implications for the energy budget 
(abstract), Bull. Am. Astron. Soc . . 22 L, 1060, 1990. 

Lindner, B. L., T. P. Ackerman, and J. B. Pollack, An efficient 
and accurate technique to compute the absorption, emission, 
and transmission of radiation by the martian atmosphere, in 

Scientific results of the NASA- sponsored s fudy . pro j eqf qn 

Mars: Evolution of volcanism. tectonics and volatiles , 

S.C. Solomon, V.L. Sharpton, and J.R. Zimbelman, ed.s, p. 
198-200, LPI Tech. Rpt. 90-06, Lunar and Planetary 

V 

Institute, Houston, 322 pp . , 1990c. 

Lindner, B. L. , T. P. Ackerman, J. B. Pollack, 0. B. Toon, and 
G. E. Thomas, Solar and IR radiation near the martian 
surface: A parameterization for CO 2 transmittance, In 

Lunar and Planetary Science XXI ■ pp . 696-697, Lunar and 

Planetary Institute, Houston, 1990. 

Lindner, B. L., J. B. Pollack, and T. P. Ackerman, C0 2 

transmittance in the Mars atmosphere: An exponential - sum 

fit for use in multiple - scattering models (abstract), 
Annales Geophvs icae . EGS special issue, p. 328, 1990. 
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CLOUDS AND CHEMISTRY ON MARS 


Bernhard Lee Lindner 


1. INTRODUCTION 

Martian O 3 shows strong seasonal and lat- 
itudinal variation, with a mixing ratio ranging 
from 2 x 10 " ppm at equatorial latitudes to a maxi- 
mum of 0 . 6 ppra over the winter polar atmosphere 
(Barth et al., 1973). The wealth of data provid- 
ed by Mariner 9 spawned many theoretical models 
to understand martian photochemistry. However, 
the polar hood (a layer of clouds at winter polar 
latitudes) has been neglected by both observa- 
tional and theoretical studies of martian aerono- 
my, even though there have been no studies that 
have shown whether these effects are indeed neg- 
ligible. Indeed, the abundance of O 3 within the 
polar hood Is virtually unknown because of the 
Inability of reflectance spectroscopy observa- 
tions of 0 3 to pierce the large optical depths 
(r) of aerosols. This study examines the varia- 
tion in the O 3 abundance on Mars with changes In 
cloud opacities from a theoretical basis; it also 
questions the efficacy of reflectance spec- 
troscopy for observing abundances on Mars. 

2. MODEL 

The O 3 abundance Is obtained by solving 
the coupled one -dimensional continuity equations 
for 0 H0 X . 0 3 . 0(*P) f 0( X D), H, OH, HO., and 
H 2°2 (discussed In detail in Lindner and Jakosky 
1985; Lindner, 1988). H 2 0 2 is allowed to con- 

dense. The surface is assumed covered with old 
ice with a albedo of 0,5. The vertical transport 
term in the one- dimensional continuity equation 
uses the eddy diffusion coefficient, K, to de- 
scribe the processes which mix the atmosphere 
vertically. Heterogeneous catalytic decomposi- 
tion of O 3 at the surface is accounted for using 
the parameterization of Kong and McElroy (1977), 
modified to account for the lower temperatures 
and snow cover at the winter poles. The radia- 
tive transfer used to calculate the intensity is 
described in detail elsewhere (Lindner, 1990; 
Lindner et al. , 1990). The discrete ordinate 
method of Stamnes et al . (1988) is used to treat 
the scattering and absorption of radiation. The 
flux is calculated for 57 solar wavelength inter- 
vals, including Rayleigh scattering and the 
absorption and scattering by O 3 , C0 2 , 0 2 , H 2 0, 
H 2 0 2 , H0 2 , cloud ice, and dust. Three types of 
clouds are frequently observed In the polar hood 
and are selected to test the interaction between 
O 3 and the polar hood; a low-lying cloud or fog, 
a cloud layer, and a hazy cloud. The vertical 
optical depth at the surface (r v ) used here is 
unity, which is typical of the polar hood. The 
single-scattering albedo Ice is taken to be unity 
and the Ucnyey-Crecns te in phase function is used 
3. OZONE AND THE POIAR HOOD 

The effect of a hazy cloud on the number 
densities of O x family members is shown in Figure 
I. H 2 0 and 1 U 0 2 ore controlled by temperature in 
the winter polar atmosphere and are not shown. 

0 3 number densities decrease at all altitudes* 
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primarily because the dominant O x loss occurs at 
onal and lat- *n altitude of 20km through H0 X , and HO actually 

tio ranging increases at this altitude. Unlike the X other 

es to a maxi- species, O 3 does not experience strong altitude- 

atmosphere dependent changes because O 3 is nearly in diffus- 

data provid- tve equilibrium (0 X Is nearly In diffusive equi- 

ical models librium, and O 3 is the predominant 0 X species 

. However, below 30km). Decreases in O 3 number densities 

winter polar naaxr the surface are nearly balanced by increases 

1 observa- in O 3 number densities at high altitudes, 

trtian aerono- H, OH, H0 2 , 0( J D), and 0( 5 P) number den- 

tudies that sities all exhibit increases above and within the 

Indeed neg- upper part of the cloud, and decreases in the 

j within the lower part of the cloud. As H, OH, 0( S P), and 

lse of the 0( X D) all have very short lifetimes, they are 

f observa- more dependent on the local photodissociation 

'al depths rates (J) than are O 3 and H0 2 , and exhibit the 

s the varia- same altitude dependence as the J's. J's In- 

\ changes In crease as much as 10 % above the cloud and within 

is Is; it also upper 10km of the cloud because of the in- 

! spec- crease In the effective solar flux resulting from 

on Mars. the scattering by the cloud. J's decrease over 

20 % in the lower part of the cloud because of a 

by solving significant reduction in the solar radiation re- 

y equations suiting from the scattering to space of the solar 

H0 2 , and flux by the upper portion of the cloud, 

and Jakosky, Changing the optical depth of the cloud 

d to con- gives results similar to those shown, except that 

d with old the magnitudes varied slightly. A cloud with r 

al transport - 3 yields virtually the same results shown, in^ 

y equation dicating that large clouds (r v - 1 - 3 ) affect mar- 

K, to de- tian aeronomy similarly and irrespective of r 

mosphere Indeed, O 3 number densities are found to vary V by 

decomposi- a t most a few percent for a cloud of any optical 

d for using depth, which agrees with Mariner 9 observations 

roy (1977), (Barth and Dick, 1974). 

peratures The most 

rhe radia- surprising result 40« , r-| 

ntenslty is Is that J's are | 

r, 1990; decreased by only f 

ordinate 20 % near the sur- \ 
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a cloud with r v - 
1 should cause far 
greater reductions 
in the available 
solar flux. The 
explanation lies 
In the large solar 
zenith angle. 
Without a cloud, 
the solar flux 
must traverse an 
effective optical 
depth of t v /p o , 
where r Is the 
vertical optical 
depth measured 
from the surface 
to infinity and /i 
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Figure 1 . Ratio of the 
number densities for the 
hazy cloud case to those 
of the cloudless case. 




Is the cosine of the solar zenith angle. In the 
winter polar regions, /i Q is quite large. If a 
cloud Is in the path of the solar rays, then most 
rays are scattered in all directions, although 
with strong forward scattering. Some solar rays 
are scattered to space and lost. However, a sig- 
nificant percentage of the solar photons are 
scattered downward at some angle 6 which is usu- 
ally less than the solar zenith angle. As a re- 
sult, r v /cos0 Is actually less than the optical 
depth traversed had the cloud not been there. 

When solar zenith angles are as large as they are 
In the winter polar regions, this is an important 
effect, explaining why the reduction in J's is 
only 20 % when a factor of 2 or so could be ex- 
pected. This effect is even more dramatic far- 
ther poleward where J's are actually larger at 
the surface when a cloud is present than they 
would be if no cloud were present. 

The overall pattern in number densities 
with a low-lying cloud or fog is similar to the 
case of the hazy cloud, only shifted toward the 
surface. Again, CU is reduced 3-4% at all 
heights. When a cloud layer is present, the J's 
and number densities develop a discontinuity at 
the altitude of the cloud layer. All short-lived 
species also show a small increase near the sur- 
face, caused by a similar increase in the J's 
near the surface. Below the cloud, the effect of 
the cloud Is to increase J's as the surface is 
approached because the solar flux passes through 
less atmosphere than it would If it traveled 
along the solar zenith angle, as explained be- 
fore. The effect is stronger for the optically 
thick wavelengths that are Important for J(C0 9 ) 
and J(H 2 0) . z 

Other latitudes and seasons are also ex- 
amined, In general, O 3 varies only a few percent 
with the occurrence of clouds, sometimes increas- 
ing slightly, sometimes decreasing slightly, 

4, DISCUSSION 

Figure 2 shows the column abundance of O 3 
observed by Mariner 9 in late winter (Barth, 

1985). While a general trend In the O 3 abundance 
with latitude is apparent In Figure 2, signifi- 
cant scatter of data points can be seen, particu- 
larly at polar hood latitudes (poleward of 40* 
latitude). Furthermore, significant variability 
occurs In the O 3 measurements made over the polar 
hood, but not in the observations made over the 
polar cap when no clouds are present (Barth et 
al., 1973). As demonstrated by the modeling re- 
sults, the presence of clouds does little to 
change the O 3 abundance, or even to change the 
altitude dependence of O 3 . Variations In temper- 
ature (cold fronts) have been claimed to account 
for much of the scatter in the data points at any 
particular latitude, because the water vapor 
abundance would vary as well (see Barth et al., 
1973; Barth and Dick, 1974; Barth, 1985). Howev- 
er, water vapor is a small source of odd hydrogen 
in the winter poplar atmosphere, and may not ac- 
count for most of the variability in Figure 2. 

Masking by clouds may also account for 
some of the observed O 3 variability, because the 
nature and opacity of the clouds in the polar 
hood change dramatically In latitude and even on 
a day-to-day basis. As the maximum O 3 abundance 
resides near the surface, spacecraft must be able 
to observe through the entire cloud In order to 
actually see the total O 3 column abundance. Tf 
reflectance spectroscopy Is used, as on Mariner 
9, then the cloud and the airborne dust must be 
traversed twice; fLrsJ. by th e Incoming solar flux 


down to the surface, and then once again upon re- 
flection from the surface out to the space- 
craft. In addition, the large solar zenith an- 
gles at winter polar latitudes mean several times 
r v of the cloud and dust must be traversed. In- 
deed, part of the observed latitudinal variation 
in O 3 may be due to the inability of the space- 
craft to observe through the increasing effective 
optical depths (r v A 0 as one goes poleward. 

This work supported by NASA grant NAStf- 

4444. 

5. REFERENCES 

Barth, C.A. , Photochemistry of the atmosphere of 
Mars, T he Photochemistry of Atmospheres. Earth. 
£he Other Planets, and Comets . (Ed. J. Levine), 
476 p., Academic Press, Orlando, Fla, 1985. 
Barth, C.A. and Dick, M.L., Ozone and the polar 
hood of Mars, Icarus . 21 , 205-211, 1974. 

Barth, C.A. , C.W. Hord, A. I. Stewart, A. L. Lane, 
M.L. Dick, and G.P. Anderson, Mariner 9 ultravi- 
olet experiment: Seasonal variation of ozone on 
Mars, Science . 122, 795-796, 1973. 

Kong, T.Y. and M.B. McElroy, The global distribu- 
tion of 0 3 on Mars, Planet. Space Sci . 21, 839- 
857, 1977. 

Lindner, B.L., Ozone on Mars: The effects of 

clouds and airborne dust, Planet. Space Sci. . 

26, 125-144, 1988. - 

Lindner, B.L. , The martian polar cap: Radiative 

effects of ozone, clouds, and airborne dust, J , 
Ccophvs. Res. . £2, 1367-1379, 1990. 

Lindner, B.L. and B.M, Jakosky, Martian atmo- 
spheric photochemistry and composition during 
periods of low obliquity, J. Ceophvs . Res f . 90 . 
3435-3440, 1985. 

Lindner, B.L. , T. Ackerman, J. Pollack, O.B. 

Toon, and G.E. Thomas, Solar and IR radiation 
near the martian surface: A parameterization 

for C0 2 transmittance. In Lunar and Planetary 
Science XXI. pp. 696-697, Lunar and Planetary 
Institute, Houston, 1990. 

Stamnes, K. , S. Tsay, W. Wlscombe, and K. 
Jayaweera, Numerically stable algorithm for 
discrete -ordinate -method radiative transfer in 
multiple scattering and emitting layered media, 
App. Optics . 21, 2502, 1988, 



Fig. 3. Mariner 9 measurements or Tim O, column 
A ltUNIMNCE DURING THE NORTHERN WINTER, L t = 330-360°. 
in Tin: Northern Hemisphere (taken with permission 
EROM Barth. 1985). 




THE EFFECT OF POLAR CLOUDS AND DUST ON THE RADIATIVE BUDGET OF THE MARTIAN POLAR CAP 


Bernhard Lee Lindner 


Atmospheric and Environmental Research, Incorporated 
Cambridge, Massachusetts. - - — ' 


1 . INTRODUCTION 

One of the most puzzling mysteries about 
the planet Mars is the hemispherical asymmetry in 
the polar caps. Every spring the seasonal polar 
cap of CO 2 recedes until the end of summer, when 
only a small part, the residual polar cap, re- 
mains. During the year that Viking observed 
Mars, the residual polar cap was composed of 
water ice in the northern hemisphere [Kieffer et 
al., 1976] but was primarily carbon dioxide Ice 
in the southern hemisphere [Kieffer, 1979]. Sci- 
entists have sought to explain this asymmetry by 
modeling Viking lander observations of atmospher- 
ic pressure (since the seasonal polar caps are 
primarily frozen atmosphere, they are directly 
related to changes In atmospheric mass). These 
models are hereafter called "polar cap/atrao spher- 
ic pressure" models because they examine the en- 
ergy balance at the surface at all latitudes and 
seasons, thus simulating the condensation and 
evaporation of CO 2 Ice over the winter polar re- 
gion. They can therefore be compared to observa- 
tions of atmospheric pressure and polar cap re- 
cession, which are. reproduced fairly well, except 
for the asymmetry In the residual polar caps 
(Leighton and Murray, 1966; Cross, 1971; Briggs, 
1974; Davies et al., 1977; James and North, 1982; 
Lindner, 1985, 1986). This paper will focus on 
how ozone, clouds, and airborne dust affect CO 2 
Ice formation and sublimation to see if they help 
explain the asymmetry in the residual polar caps. 

2. MODELING PROCEDURE 

The radiative flux striking the surface 
Is calculated for 57 solar and 10 infrared wave- 
length intervals from 0 to 100 fim including the 
absorption, scattering, and emission by 0^, CO 2 , 
clouds, and dust (see Lindner, 1990). Ultravio- 
let absorption by O 2 , fyO, HO 2 , and H 2 O 2 Is also 
Included, as Is Rayleigh scattering. The dis- 
crete ordinate method of Stamnes et al. (1988) 
treats the scattering, emission, and absorption 
of monochromatic radiation. The exponential sum 
method allows CO 2 transmittance to be Incorpora- 
ted in a scattering model [see Lindner et al., 
1990]. The albedo for the polar cap is 0.5 at 
solar wavelengths [Kieffer, 1979; James and 
Lumme , 1982]. The IR albedo of the polar cap is 
assumed to be zero [Kieffer, 1970; Sraythe, 1975; 
Wiscorabe and Warren, 1980] . Three latitude cases 
were studied: 57*N, 70*N, and 90*N, and late win- 
ter conditions were assumed (L s - 343*). These 
three cases effectively simulate the edge of the 
polar cap, the edge of polar night, and polar 
night, respectively, for all solar longitudes and 
both polar caps. The temperature profile at 57*N 
latitude rises linearly with altitude from 150 K 
at the surface to 170 K at 6 km altitude and then 
falls to 130 K at 40 km [Lindal et al. , 1979]. 

The temperature profile at 70*N had a peak tem- 
perature of 155 K at 6 km, and polar night tem- 
peratures fell linearly from 150 K at the surface 


to 130 K at 40 km. The value used for albedo and 
temperature had little effect on our conclusions. 

3. THE EFFECT OF CLOUDS AND DU'ST 

Sublimation depends on the total flux 
that ls absorbed, which is a function of the cap 
albedo. Weighting the fluxes we calculate with 
the coalbedo results In the absorbed flux, as 
shown In Table 1. Almost no change in the ab- 
sorbed flux occurs over the range that dust and 
cloud optical depth experiences in the atmo- 
sphere, a result also found by Davies [1979] and 
Paige [1985]. The increase In the infrared flux 
due to thermal emission by dust balances the de- 
crease in the solar flux due to dust absorption. 

It Is at the polar night latitudes that 
the radiative effects of clouds and background 
dust are most significant, as shown in Table 2. 

IR fluxes absorbed by the polar cap increase up 
to an order of magnitude as the cloud opacity and 
dust opacity increase. Even simply a cloud opti- 
cal depth of 0.5 and dust optical depth of 0.2 
would result in the sublimation of the order of 
40 cm of CO 2 ice over 1 martian year. When in- 
serted in a model, this would accelerate the pre- 
dl cted loss of al l the CO 2 ice on the northern 


TABLE 1. Flux Absorbed by the Polar Cap Near 
Its Edge for Various Cloud and Dust 
Opacities ( 5 7 * N Latitude, L s - 343*) 
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0.0 

0.2 

0.5 
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polar cap by approximately 40 sols (based on Ice 
depth vs. season from Leighton and Murray [1966], 
Cross [1971], Briggs [1974], Davies et al [1977], 
Lindner [1985], and Jakosky and Haberle [1990]), 
While polar cap/atmospheric pressure mo- 
dels have been unable to explain the dichotomy of 
the residual polar caps, these models have been 
fairly successful at reproducing observations of 
atmospheric pressure and polar cap recession up 
to 75 or 80* latitude. The radiative effects of 
clouds and background dust have the added benefit 
of not significantly affecting model predictions 
of the annual variation in atmospheric pres- 
sure. Clouds and background dust have a neutral 
effect along the cap edge, becoming more impor- 
tant as one approaches the pole, particularly 
within 80* latitude . The integral of seasonal 
COp ice over the 80*S to 90*S latitude band 
(where the residual polar cap exists) is an order 
of magnitude less than the integral of seasonal 
COj ice over the planet. Hence, a process which 
affects primarily the residual polar cap area 
will not significantly affect predictions of at- 
mospheric pressure. Also, the rate of recession 
of the seasonal polar cap would remain essential- 
ly the same as that predicted by polar cap/atmo- 
spheric pressure models equatorward of 75* lati- 
tude, where good agreement with observations was 
previously obtained. The behavior near the pole 
would be significantly modified, making for bet- 
ter agreement with observations than previously 
obtained. The rate of recession of the polar cap 
is an important test, especially since it repeats 
reliably over several years of observation [James 
and Lurame, 1982; James, 1982]. 

4. IMPLICATIONS FOR THE POUR CAP ASYMMETRY 

The radiative effects of clouds and back- 
ground dust would be most noticeable near the re- 
sidual polar cap, which spends all of the winter 
in polar night. Observational evidence suggests 
a greater abundance of clouds over the northern 
polar cap than over the southern polar cap. 
Moreover, it Is also expected that the formation 
and sublimation phases of the seasonal southern 
polar cap will occur with higher background dust 
opacities than the formation and sublimation 
phases of the seasonal northern polar cap. Both 
of these asymmetries would preferentially in- 
crease the fluxes absorbed by the residual polar 
cap in the north. 

Furthermore, the surface pressure is lo- 
wer over the southern residual polar cap than the 
northern residual polar cap. If mixing ratios of 
cloud and dust are assumed to be the same for 
both poles, then the pressure differential be- 
tween the poles would result In a dust and cloud 
opacity over the residual polar cap in the south 
that is 30% less than over the north during the 
crucial winter months. Mixing ratios of dust may 
actually be lower' over the southern polar cap 
during winter months because the mechanisms for 
raising and maintaining atmospheric dust are 
pressure dependent [Pollack and Toon, 1982]. 

Lower mixing ratios over the southern pole would 
accentuate the asymmetry in cloud and dust opaci- 
ty. These hemispherical asymmetries would allow 
relatively less frost to accumulate in the north 
during the fall and winter and cause the CO 2 sur- 
face Ice to sublimate more rapidly in the north 
during the following spring and summer. 

However, the effects of clouds and back- 
ground dust would not allow CC>2 ice to survive 


year-round in the south, as is observed. Other 
processes which might allow that include penetra- 
tion of radiation into and through the ice bidi- 
rectional surface reflectance, ice albedo depen- 
dence on solar zenith angle, season, and latitude 
[Vis combe and Warren, 1980; Squyres and Veverka 
1982; Paige 1985], lower ice emisslvity (S. War- 
ren, personal communication, 1988), decreased 
heat conduction from the surface due to existing 
Ice [Jakosky and Haberle, 1990], snowfall [Pol- 
lack and Haberle, 1988], surface residues [Saund- 
ers et al., 1986], surface roughness, and wind 
shifting of ice [Briggs, 1974; James et al., 

1979; Kleffer, 1979], The importance of each of 
these processes Is currently being assessed by 
B.L. Lindner and B.M. Jakosky. Most of these 
processes would lengthen CO 2 Ice survivability 
either by increasing condensation or decreasing 
sublimation of CO 2 Ice and would counter the in- 
creased sublimation caused by clouds and dust. 
However, most of these processes work equally ef- 
fectively at both poles and do not generate an 
asymmetry. Hemispherical asymmetries in cloud 
and dust opacity and in surface elevation would 
result In a hemispherical asymmetry in the net 
radiative flux from clouds and dust which is ab- 
sorbed by the polar cap, helping explain the 
asymmetry in the polar caps. 

This work paid by NASA contract NASW-4444. 
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Abstract . The solar and thermal flux striking 
the polar cap of Mars is computed for various 
ozone , dust, and cloud abundances and for three 
solar zenith angles. Ozone does not significant- 
ly affect the total energy budget of the polar 
cap. Hence the observed hemispherical asymmetry 
in ozone abundance causes only an insignificant 
hemispherical asymmetry in the polar caps. Ver- 
tical optical depths of dust and cloud ranging 
from zero to 1 cause little change in the total 
flux absorbed by the polar cap near Its edge but 
increase the absorbed flux significantly as one 
travels poleward. Hemispherical asymmetries in 
dust abundance, cloud cover, and surface pressure 
combine to cause a significant hemispherical 
asymmetry in the total flux absorbed by the re- 
sidual polar caps, which helps to explain the 
dichotomy In the residual polar caps on Kars. 
Other processes which affect the energy budget of 
the polar cap are proposed and reviewed, particu- 
larly with respect to their Interaction with the 
radiative effects of clouds and dust. 

Introduction 

One of the most puzzling mysteries about the 
planet Kars is the hemispherical asymmetry In the 
polar caps. Every spring the seasonal polar cap 
of <X >2 recedes until the end of summer, when only 
a small part, the residual polar cap, remains. 
During the year that Viking observed Mars, the 
residual polar cap was composed of water ice in 
the northern hemisphere [Kieffer et al., 1976] 
but was' primarily carbon dioxide Ice In the 
southern hemisphere [Kieffer, 1979]. Scientists 
have sought to explain this asymmetry by modeling 
Viking lander observations of atmospheric pres- 
sure (since the seasonal polar caps are primarily 
frozen atmosphere, they are directly related to 
changes In atmospheric mass). These models are 
hereafter called "polar cap/atmospheric pressure" 
models because they examine the energy balance at 
the surface at all latitudes and seasons, thus 
simulating the condensation and evaporation of 
C (>2 ice over the winter polar region, and can 
therefore be compared to observations of atmo- 
spheric pressure and polar cap recession [Leigh- 
ton and Murray, 1966; Cross, 1971; Briggs, 1974; 
Davies et al., 1977; James and North, 1982; Lind- 
ner, 1985, 1986]. Polar cap/atmospherlc pressure 
models reproduce most aspects of the observed an- 
nual variation in atmospheric pressure fairly ac- 
curately. Furthermore, the predicted recession 
of the northern polar cap in the spring agrees 
well with observations, including the fact that 
the CO 2 ice is predicted to completely sublime 
away. However, these models all predict that the 
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carbon dioxide Ice will also sublLme away during 
the summer in the southern hemisphere, as shown 
in Figure 1. This paper will focus on how ozone, 
clouds, and airborne dust affect CO^ Ice forma- 
tion and sublimation to see if they help explain 
the hemispherical asymmetry In the residual polar 
caps . 

Observations of Ozone . Clouds . and A irborne — Pu$£ 

Ozone has been observed on Mars by Mariner 6 
and Mariner 7 [Barth and Hord, 1971; Lane et al., 
1973], by Mariner 9 [Barth et al., 1973; Lane et 
al., 1973; Barth and Dick, 1974], by MARS 5 
[Krasnopol # skll et al. , 1975], and by Earth-based 
experiments [Broadfoot end Wallace, 1970; Noxon 
et al., 1976; Traub et al.,‘l979] v Global-aver- 
age ozone abundances Integrated from the surface 
to space are of the order of 1 pm atm (1 micron 
atmosphere w 1 pm atm — 2.69 x 10 cm ). Ozone 
experiences strong latitudinal variation. In- 
creasing significantly over the winter polar 
regions. In the northern hemisphere, the maximum 
abundance of ozone observed was 60 pm atm near 
60 # N latitude in late winter, while the maximum 
abundance observed In the southern hemisphere was 
30 pm atm. This hemispherical asymmetry in ozone 
is reproduced by theoretical modeling [Kong and 
McElroy , 1977; Shimazaki and Shimizu, 1979; Lind- 
ner, 1985]. There is also strong day-to-day var- 
iability In ozone, which arises from the variabi- 
lity in temperature and In the polar hood [Barth 
and Dick, 1974; Lindner, 1988]. 

The polar hood, a blanket of clouds which 
covers the winter polar region, has been observed 
on a latitudinal and seasonal basis by Mariner 9, 
MARS 3, and Viking [Leovy et al., 1972; Masursky 
et al., 1972; Briggs and Leovy, 1974; Martin, 
1975; Moroz, 1976; Briggs et al., 1977; Anderson 
and Leovy, 1978; Tillman et al. , 1979; Kondratiev 
and Hunt, 1982; Kahn, 1984; Christensen and 
Zurek, 1984]. The following synopsis is extrac- 
ted from these papers. The polar hood Is first 
observed at high latitudes in the northern hemi- 
sphere at about the time of autumnal equinox, 
gradually spreading to obscure the surface to 
latitudes as low as 40*, and then retreating to 
the pole in the spring. The disappearance of the 
hood occurs at about the vernal equinox. At 

lower latitudes, the polar hood has strong daily 

variation and highly variable and often complex 
vertical structure and is diffuse on scales of 
less than 1 km [Briggs and Leovy. 1974]. The 

main bulk of the polar hood Is below 20 km alti- 

tude. The thickness at one point is estimated at 
1-2 km. centered at 15 km height, with an average 
mixing ratio of ice of the order of 10 g 
cm’ . The hood is often a wave cloud pattern at 
lower latitudes, and clouds 10 km high with half 
an optical depth can appear at the hood's edge. 
Equatorward of 60* latitude in winter the compo- 
sition of the clouds may be primarily water 
ice. Poleward of that, the hood is more diffuse. 
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Fig. 1. The seasonal recession of the south 
polar cap as observed over the last 20 years 
[James and Lumme , 1982 J and as predicted by Leigh- 
ton and Murray (1966), Briggs [1974], Davies et 
al. [1977}, Y. Naruml, (unpublished data, 1980) 
and James and North [1982}. (The aerocentric lon- 
gitude of the Sun, L s , Is the seasonal Index; L fi — 
0*^, 90*, 180*, and 270* correspond to northern 
spring equinox, summer solstice, autumnal equinox, 
and vlnter solstice, respectively.) 

lower to the ground, and may be composed of car- 
bon dioxide ice. In general, the polar hood ls 
diffuse and structureless on scales smaller than 
a few kilometers. 

Observers do note a reduction In polar hood 
clouds during global dust storms, although dust 
obscuration of most of the lower atmosphere 
causes the degree of any reduction In cloud cover 
to be Inconclusive. Imaging observations of the 
polar hood in the polar night are nonexistent. 
However, COj ice clouds have been Inferred to 
exist in the polar night from Viking Infrared 
Thermal Mapper (IRTM) data [Hunt et al., 1980; 
Paige, 1985}. Morning fog Is observed regularly 
over the Viking landers, and observations for 
several seasons have been reported [Pollack et 
al., 1977; Ryan and Sharman, 1981; Ryan et al., 
1982; Jakosky, 1985}. However, morning fog is 
not expected at winter polar latitudes, since 
winter polar latitudes experience little diurnal 
temperature change. 

Clouds are observed to be more prevalent over 
the edge of the northern polar cap than over the 
edge of the southern polar cap, although coverage 
is far from complete [James, 1983; Christensen 
and Zurek, 1983; James et al., 1987]. Clouds 
could not be seen in polar night with the imagers 
sent to Mars to date, although their existence 
has been suggested based on analysis of Viking 
Infrared data [Hunt et al., 1980; Paige, 1985]. 
Observations in polar night are not able to dis- 
cern whether or not any of the hemispherical 
asymmetry that may exist along the edge of the 
polar cap exists in polar night as well, although 


there is no evidence to suggest any other behav- 
ior . 

The optical depth of dust in the atmosphere 
also experiences significant seasonal varia- 
tion. Global dust storms tend to occur during 
southern spring, when Mars is near perihelion 
[Martin, 1984; Tillman and Leovy, 1984} , and tend 
to have vertical optical depths of dust of 1.0 or 
greater at visible wavelengths [Pollack et al., 
1979; Thorpe, 1981]. During the remainder of the 
year, atmospheric dust opacities are substantial- 
ly lower, although still appreciable, and are 
hereafter called background dust. The vertical 
optical depth of background dust was observed to 
vary from 0.3 to 1.0 In the visible as determined 
by imaging the Sun from Viking landers 1 and 2 
[Pollack et al., 1979]. Similar background dust 
optical depths were inferred from contrast obser- 
vations by Viking orbiter imaging [Thorpe, 1981], 
Inferred from Viking lander pressure oscillations 
[Zurek, 1981] , deduced by studying polar cap 
reflection from Earth-based observations [Lumme 
and James, 1984], and Inferred from modeling the 
annual variation In pressure [James and North, 
1982} . 

The frequency and Intensity of dust storms 
have apparently varied greatly from year to year 
over the past 60 years [Martin, 1984]. The se- 
cond and third Kars years that Viking observed 
were relatively calm, while the fourth year had 
the most intense global dust storm observed by 
Viking, during which Viking lander 1 failed 
[Tillman and Leovy, 1984], Earlier Mariner 9 
imaging of surface contrast noted significantly 
less dust than observed by Viking, with optical 
depths of dust ranging from 0.1 to 2 [Leovy et 
al., 1972; Masursky et al., 1972]. The same var- 
iation in dust optical depth was Inferred from 
Mariner 9 infrared interferometer spectrometer 
(IRIS) spectra, at 3000 A [Toon et al., 1977] and 
was determined by the reduction In the solar flux 
reflected by the polar cap, as observed by the 
Mariner 9 UV spectrometer [Pang and Hord, 1973]. 

Observations are Inadequate to accurately 
describe any latitudinal variation in background 
dust opacity; orbiter data fail to note any sig- 
nificant latitudinal variation, other than a con- 
sistently lower dust opacity over both winter 
poles. Dust opacities In the winter polar atmo- 
sphere have been observed to be a factor of 2 or 
so lower than elsewhere over the planet when the 
global dust storms exist [Masursky et al., 1972; 
Leovy et al., 1972; Pang and Hord, 1973; Kieffer, 
1979; James et al., 1979; Lumme and James, 1984; 
Jakosky and Martin, 1987], and the same may be 
true during background dust conditions. Surface 
dust Is covered by ice at winter polar latitudes, 
hindering that source of dust; the winter polar 
atmosphere is stable, facilitating dust settling; 
and C (>2 snow forms using dust particles as nu- 
clei, further accelerating dust removal [Pollack 
et al., 1979], All of these processes may be 
active during background dust conditions, result- 
ing in dust opacities that may be becoming gradu- 
ally lower as one approaches the winter pole. 

However, the dust seasonal cycle results in a 
background dust opacity over the Viking lander 
which is higher from L s - 180* to 360* (southern 
winter and spring) than from L s - 0* to 180* 
(northern winter and spring) [Pollack et al., 
1979]. Assuming that little hemlspber lea 1 asyrn- 
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raetry occurs for the same time frame in back- 
ground dust opacity [Martin, 1986], then the for- 
mation and sublimation phases of the southern 
polar cap will occur with higher background dust 
opacities than those of the northern polar cap. 
As a result, background dust could affect the 
formation and sublimation of the southern cap 
more than the northern cap at polar night lati- 
tudes . 

Previous Work in This Area 

Kuhn et al. [1979] suggested that ozone could 
play an Important role in the energy budget of 
the polar cap. Kuhn et al. showed that ozone 
heating of the winter polar atmosphere was non- 
negllgible compared to COj heating, particularly 
near 10 km altitude, and could therefore affect 
any CO2 condensation in the atmosphere. However, 
atmospheric heating by aerosols in the winter 
polar atmosphere has since been shown to be 
greater than either or CO2 heating, particu- 
larly near 10 km altitude [Lindner and. Thomas, 
1983; Lindner, 1985] . Thus, ozone does not make 
a large enough contribution to the atmospheric 
temperature to significantly affect downwelling 
thermal radiation which strikes the polar cap or 
any condensation of CO 2 in the atmosphere. Ozone 
also affects the surface energy budget by absorb- 
ing incoming solar radiation which would other- 
wise strike the surface. This issue is addressed 
in the results section of this paper. 

Clouds will also affect the energy budget by 
scattering solar radiation and emitting thermal 
radiation. The radiative effects of clouds were 
shown to be important by Briggs [1974] and James 
and North [1982]. However, the magnitude of this 
importance was not accurately determined, since 
they ignored the effects of clouds at solar wave- 
lengths. Furthermore, they included infrared ef- 
fects with an atmospheric emissivlty of 0.3 when 
clouds existed, an emissivlty of 0.15 for cloud- 
less conditions, and a single atmospheric temper- 
ature, assumptions which ignore the variations in 
cloud opacity and temperature which occur. 

The attenuation of solar radiation due to 
absorption and scattering by dust and the in- 
crease in thermal radiation due to emission by 
dust will also influence the energy budget of 
polar cap, as previously suggested by Briggs 
[1974], Davies [1979], James et al. [1979], Kief- 
fer [1979], Hartin and Kieffer [1979], James and 
North [1982], and Paige [1985]. These studies 
all examined the radiative effects of the global 
dust storms, which occur primarily in southern 
spring, and showed that global dust storms do not 
appear to significantly affect polar cap reces- 
sion. However, appreciable amounts of background 
dust exist year-round. Davies [1979], James and 
North [1982] , and Paige [1985] examined the ra- 
diative effects of background dust and noted a 
neutral contribution on the whole (a slight de- 
crease in some cases and a slight increase In 
others) . 

However, these background dust studies all 
ignored polar night latitudes (those latitudes 
and seasons where the Sun remains below the hori- 
zon ail day). In addLtion, Paige [1985] studied 
only the 90" latitude In his study. Furthermore, 
all studies ignored the overlap of the radiative 
effects of CO ^ and dust In the 15-pra wavelength 


region, where both CO2 and dust emit the most 
thermal radiation [Lindner, 1985], as well as the 
overlap of CO2 and dust at near- inf rared (NIR) 
wavelengths. Improper treatment of this overlap 
will result In errors in the downwelling thermal 
radiance of the order of 30% [Lindner, 1985]. 
Moreover, Davies [1979] and James and North 
[1982] accurately modeled the radiative effects 
at solar wavelengths, but Included thermal heat- 
ing of the surface by simply assuming that half 
of all the solar radiation absorbed by dust is 
radiated thermally to be absorbed by the sur- 
face. This hypothesis assumes that dust is the 
primary source of radiative heating and cooling; 
however, COj heating and cooling is also appre- 
ciable during background dust conditions [e.g., 
Lindner, 1985]. Furthermore, atmospheric absorp- 
tion of planetary thermal emission is also appre- 
ciable. A better approach to including IR heat- 
ing of the polar cap is to compute the thermal 
emission based on the observed atmospheric tem- 
peratures, €<>2 opacities, and dust opacities, as 
Is done here. 

Modeling Procedure. 

The radiative flux striking the surface is 
calculated for 57 solar and 10 infrared wave- 
length Intervals from 0 to 100 pm including the 
absorption, scattering, and emission by 0^, CO2 , 
clouds, and dust. Ultraviolet absorption of 
solar flux by (>2. H^O, HO2, ®2^2 a ^*° 
included, as is Rayleigh scattering. Atmospheric 
composition is taken as 95% CO2 and 0.13% O2 
[Owen et al., 1977]. Season-dependent C0 2 abun- 
dances are taken from Hess et al. [1980], after 
correcting for elevation [Llndal et al. , 1979; 
Jakosky and Farmer, 1982] and accounting for pos- 
sible circulation- induced pressure gradients 
[Haberle et al., 1979]. Ozone abundances were 
taken from Barth et al. [1973], assuming that 
ozone mixing ratios in the polar night were the 
same as at the edge of polar night and adopting 
altitude profiles from Lindner [1988]. H2O 
abundances are based on the work of Jakosky and 
Farmer [1982] and Lindner [1988], and HO^ and 
H2O2 abundances are taken from Lindner [ 1988 J . 

Transmission functions for the 2.7-, 4.3-, and 
15.0 -pm bands of carbon dioxide are taken from 
the line-by-line model results of Gal'tsev and 
Osipov [1979]. The transmission function, Tr, as 
a function of temperature T, pressure P, and CO2 
column abundance U was extrapolated from the 
Gal'tsev and Osipov results (subscript C) to tem- 
peratures below 200 K by 

Tr(T.P,D) - 1-[(1-Tr c (200 K.P.U)) (T/200 K) Q ] 

The exponential Q was found to be 0.45, 0.3, and 
0.8 for the 2.7-, 4.3-, and 15 -pm bands, respec- 
tively, when the temperature dependencies for the 
Pollack et al. [1981] transmission functions were 
recast in this form. Using a modified version of 
the FASCOD transmission model (Clough et al. , 
1986), the accuracy of these transmission func- 
tions was confirmed, and transmission functions 
were obtained for the 1.3-, 1 . 4 - , 1.6- % 2.0-, 
4.8-, and 5.2-pm bands of CO2 (B. L. Lln<iner et 
al., manuscript in preparation, 1989). 

The properties of clouds and dust are extrapo- 
lated from observations at mid- lat itudes , given 
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the lack of observations at winter polar lati- 
tudes. Background dust opacities over winter 
polar latitudes may be less than the Viking 
lander latitudes, where they vary from 0.2 to 1.0 
[Pollack et al., 1979; Lumrae and James, 198A] . 
However, vertical optical depths of background 
dust vary from 0 to 1.0 in this work to account 
for all possible scenarios. The wavelength de- 
pendence of the dust opacity is taken from Toon 
et al. [1977]. A Caussian profile describes the 
vertical distribution of dust, with dust opaci- 
ties confined mostly below 20 km altitude for 
normal conditions and 50 km altitude for dust 
storm conditions [see Anderson and Leovy, 1978; 
Zurek, 1982]. The single-scattering albedo of 
airborne dust as a function of wavelength is 
taken from Zurek [1978, 1982] and Toon et al. 

[1977] at solar and infrared wavelengths, respec- 
tively, using a solar average of 0.86 [Pollack et 
al., 1979], The Haze-L phase function is used to 
describe the scattering of radiation by Martian 
dust [Toon et al., 1977]. The emisslvity of air- 
borne dust is high and has been calculated as a 
function of wavelength from theory and observa- 
tions [Toon et al., 1977; Simpson et al., 1981]. 

Clouds were assumed to be diffuse, with a dis- 
tribution in altitude of the mixing ratio of 
cloud ice like that assumed for airborne dust 
[Briggs and leovy , 1974], Cloud opacities are 

observed to vary significantly, and a range of 0 
to 3 Is used in this work. The a ingle -scattering 
albedos of both C0 2 and HjO ice clouds are high 
and are assumed to be unity at all wavelengths 
from DV to IR. This agrees with single -scatter- 
ing albedos for Martian water ice clouds derived 
from IRTM data [Hunt, 1979] and derived for CO2 
ice clouds from Mie theory [Hunt et al. , 1980J 

^ or *11 wavelengths, to the extent that ice par- 
ticle radii are known. Single-scattering albedos 
for the clouds are also given by Moroz [1976] and 
Freeman and Llou [1979]. These observations 
agree well with theoretical calculations [WIs- 
combe and Warren. 1980]. The analytic phase 
function which has been used most extensively in 
the literature was introduced by Henyey and 
Greenstein [19A1] . The Henyey-Greenstein phase 
function has been shown to be valid for ice [Han- 
sen, 1969], and asymmetry factors as a function 
of wavelength from the UV to IR have been derived 
for Martian clouds, with a great deal of uncer- 
tainty [Hunt, 1979; Hunt et al., 1980]. 

The discrete ordinate method of Stamnes and 
Conklin [198A] treats the scattering, emission, 
and absorption of monochromatic radiation through 
the Martian atmosphere. The exponential sum 
method allows the banded wavelength structure of 
C0 2 in the infrared (IR) and near- Infrared (NIR) 
to be treated as monochromatic, allowing for its 
incorporation In a scattering model [see Freeman 
and Liou, 1979J. Eight terms in the exponential 
sum were found to give a good approximation to 
the transmission In the 15-/im band of C0 2 , while 
four terms were used for the NIR bands of C0 2 (B. 

L. Lindner et al.. manuscript In preparation, 
1989). The use of the exponential sum technique 
allows the overlap of C0 2 and dust opacities in 
the 15 -pm and NIR wavelength regions to be pro- 
perly treated. 

Atmospheric properties are zonally averaged 
and assumed azimuthally independent. The region 
from the surface to A0 km altitude is broken into 


20 2-km-thick layers to account for vertical in- 
homogeneity. The Chapman function Is used to 
approximate the slant path in place of the secant 
function [e.g.. Smith and Smith, 1972), because 
the winter polar atmosphere always has large 
solar zenith angles, and the secant function Is 
In error for large angles. The model Is diumal- 
ly averaged [e.g., Cogley and Borucki, 1976). An 
albedo for the polar cap at solar wavelengths 
3-5.4 pm) of 0.5 is used, which is appropriate 
for late northern winter [Kieffer, 1979; James 
and Lumme, 1982]. The IR (5.A-100 pm) albedo of 
the polar cap Is assumed to be zero [Kieffer, 
1970; Smythe, 1975; WIscombe and Warren, 1980]. 
Three latitude cases were studied: 57*N, 70*N. 

and 90 N, and late winter conditions were assumed 
" 3A3*). These three cases effectively simu- 
late the edge of the polar cap, the edge of polar 
night, and polar night, respectively, for all 
solar longitudes and both polar caps. 

The temperature profile at 57*N latitude that 
is used In this study rises linearly with alti- 
tude from 150 K at the surface to 170 K at 6 km 
altitude and then falls with increasing altitude 
to 130 K at A0 km. This temperature profile is 
extracted from radio occulation observations made 
near winter polar latitudes during background 
dust conditions [Lindal et al., 1979]. These 
temperatures also agree with those obtained by 
the Viking IRTM experiment [Kieffer, 1979; Mar- 
tin, 1984]. The diumally averaged heating and 
cooling rates calculated using these temperatures 
were very nearly in balance. This Indicates that 
the winter polar atmosphere Is nearly In radia- 
tive equilibrium, which Is the assumption that Is 
used for calculating atmospheric temperatures 
farther poleward. The temperature profile at 
70*N latitude was the same as that used at 57*N 
latitude except that it had a peak temperature of 
155 K at 6 km. Polar night temperatures were 
assumed to have no inversion and fell linearly 
from 150 K at the surface to 130 K at 40 km. 
Temperatures are markedly higher at the time of 
the global dust storms [e.g.. Pollack, 1978], but 
global dust storms are not considered here. 


Radiative Fluxes at the Polar Cap Surface 

Table 1 shows the radiative fluxes from 0 to 
100 n Integrated over UV, visible, NIR, and 
three IR wavelength bands. 0 3 absorbs almost two 
thirds of the available solar flux from 0 to 0.3 
!*m at 57*N latitude. (The flux in the 0- to 0.3- 
/im bin is 6 J cm sol for no dust and no 0 3 
scenarios for this latitude and season. 1 sol ■ 
1 Martian day.) However, ozone absorption has 
only a minor effect on the overall flux at the 
surface, due to the small contribution that solar 
radiation In the UV makes to the overall energy 
budget. Very little 0 3 or C0 2 absorption occurs 
in the 0.3- to 1.0 -fim wavelength bin. Absorption 
In the C0 2 NIR bands prevents approximately 5% of 
the available solar flux from 1 to 5 from 
striking the surface at 57’N latitude and L - 
343 . The 15 -/im band of C0 2 is the dominant gas- 
eous source of thermal radiation which strikes 
the polar cap [see Lindner, 1985]. The IR «*lux 
from the Sun can be seen in the 5. A- to 10-/im and 
21- to 100-/im wavelength bins for the dust-free 
case. 
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TABLE 1. Flux Striking the Surface Near the Edge of the Polar Cap (57"N Latitude, L s - 343") 


Wavelength 
Interval, fx m 


Dust 

Opacity 


0.0 

0.2 

0.5 

1.0 

0 . 0 - 0 . 3 

2.132 

1.429 

0.823 

0.364 

0 . 3-1.0 

455.608 

391.262 

321.688 

247.947 

1.0-5. 4 

194.296 

167.001 

138.951 

109.501 

5.4-10.0 

2.221 

3.255 

4.361 

5.315 

10 . 0 - 21.0 

19.894 

41.189 

61.311 

79.603 

21 . 0 - 100.0 

0.044 

16.603 

36.640 

62.910 

Total flux 

674.195 

620.739 

563.774 

505.640 

Absorbed flux 

348.177 

340.893 

333.043 

326.734 


The vertical optical depth of dust is given. 


Flux values are given In units of J cm 


sol 


<1 


Dust affects the energy budget of the polar 
cap by absorbing and scattering to space solar 
radiation that would otherwise strike the cap and 
by emitting thermal radiation , some of which 
strikes the cap. Dust absorbs and scatters 
strongly In the UV and moderately tn the visible 
and NIR as seen in Table 1. Dust also emits 
thermally at all wavelengths. Note that the flux 
which strikes the surface from thermal emission 
by even small abundances of airborne dust ls of 
the same order as thermal emission in the 15-ywn 
band of COj. Clearly, both CO 2 end dust thermal 
emission are important for the energy budget of 
the polar cap. It ls also Important to note that 
in the 15-/«n wavelength band, both dust and CO 2 
emission are important. Hence, it is important 
to treat dust and CO 2 simultaneously as ls done 
here, and not to calculate cooling rates Individ- 
ually and sum them. 

The total wavelength- Integrated flux decreases 
with increasing dust optical depth, as the in- 
crease in thermal radiation does not equal the 
decrease in solar flux. However, sublimation 
depends on the total flux that ls absorbed, which 
is a function of the cap albedo. Vhereas the 
typical albedo of old ice is 0.5 at solar wave- 
lengths [Kieffer, 1979; James and Lumme, 1982], 
the albedo at IR wavelengths is close to zero 
[Kieffer, 1970; Smythe, 1975; VIscombe and War- 
ren, 1980]. Weighting the flux with the coalbedo 
(1-albedo) results In the absorbed flux, given at 
the bottom of Table 1. Almost no change In the 
absorbed flux occurs over the range that dust 
optical depth experiences in the Mars atmosphere 
(omitting dust storms), a result also found by 
Davies [1979] and Paige [1985]. The increase in 
the infrared flux due to thermal emission by dust 
balances the decrease In the solar flux due to 
dust absorption, although not completely, result- 
ing in a slight decrease in absorbed flux. Glo- 
bal dust storms are another matter, as the atmo- 
spheric temperature and the meridional transport 
of heat are much higher. However, the effect of 
dust storms on polar cap recession was studied by 
Briggs [1974], James and North [1982], and Paige 
[1985], and dust storms do not appear to offer an 
explanation for the persistence of CO 2 ice In the 
southern hemisphere. 


As shown in Table 2, all solar fluxes at 70* N 
latitude are lower than at 57 *N -latitude due to 
the Increased solar zenith angle and decreased 
daylight hours. The larger slant path allows O 3 , 
COj , and dust to be more effective at preventing 
solar flux from reaching the surface; effective 
optical depths are nearly an order of magnitude 
larger than vertical optical depths. Thermal 
fluxes at 70*N latitude are also lower than ther- 
mal fluxes at 57 *N latitude (Table 1) because 
atmospheric temperatures are lower [Lindal et 
al. , 1979]. However, the net effect on the flux 
absorbed by the polar cap Is the same as at 57*N 
latitude. In that changes in dust opacities from 
0 to 0.5 cause virtually no change^ In absorbed 
flux. 

However, while dust has little Influence on 
the total flux absorbed by the polar cap at sun- 
lit latitudes, a noticeable result occurs in the 
polar night region (that part of the planet where 
the Sun does not Illuminate the surface at all 
during an entire day) . Table 3 shows the thermal 


TABLE 2. Flux Striking the Surface of the 
Polar Cap Near the Edge of Polar Night 
(70*N Latitude, L s - 343') 


Wavelength 
Interval, /im 


Dust Opacity 


0.0 

0.2 

0.5 

0 . 0 - 0 . 3 

1.105 

0.552 

0.264 

0 . 3-1.0 

220.890 

164.544 

121.753 

1.0-5. 4 

93.550 

73.743 

58.566 

5.4-10.0 

1.078 

1.731 

2.432 

10 . 0 - 21.0 

16.348 

32.923 

48.696 

21 . 0 - 100.0 

0.021 

14.497 

32.051 

Total flux 

332.992 

287.990 

263.762 

Absorbed flux 

175.220 

168.571 

r v 173.471 

The vertical 

optical 

depth of dust ij 

; given. 

Flux values are 

given In 

units of J cm'^ 

sol' 1 . 
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TABLE 3. Flux 
Cap in Polar 

Striking 
Night (90 

the Surface of the 
*N Latitude, L s - 

Polar 

343*) 



Dust Opacity 


Wavelength 
Interval, /am 

0.0 

0.2 

0.5 

5.4-10.0 

10.0-21.0 

21.0-100.0 

0.000 

14.714 

0.000 

0.708 

29.441 

13.511 

1.391 

43.465 

29.909 

Absorbed flux 

14.714 

43.660 

74.765 


The vertical optical depth of dust is given. 
Flux values are given in units of J cm - * sol . 


flux absorbed by the polar cap for three dust 
opacities (solar fluxes are obviously zero). 
Again, note that thermal fluxes are lover at 
polar night latitudes than at sunlit latitudes 
» Tables 1 and 2) due to the extremely cold 
temperatures. Even so, there is a substantial 
Increase in absorbed flux with dust opacity, even 
for very small opacities. For comparison, the 
thermal emission by the surface to space is only 
260 J cm ^ sol asstiming unit emissivity. In 
terms of a sublimation rate, 10 J cm - ^ sol“^ 
would sublimate approximately 0.02 cm of C0 2 ice 
P e 5- sol (this assumes an ice density Of 1 g 
Cm .) Hence, a period of 1 Martian year with a 
vertical optical depth of dust of 0.2 would re- 
sult in a substantial loss of the order of 20 cm 
of C0 2 ice. 

The Effect of Clouds 

Clouds also affect the energy budget of the 
polar cap. At 57*N latitude, the radiative ef- 
^®ots on the absorbed flux of a combination of 
clouds and dust are fairly neutral, with a slight 
downward bias (Table 4). Clouds slightly de- 
crease the flux absorbed by the polar cap for 
small dust opacities (rD - 0 or 0.2), with the 
extent of the decrease proportional to the cloud 
opacity. However, clouds actually slightly in- 
crease the flux absorbed by the polar cap at 57*N 


latitude for moderate or large dust opacities (rD 
- 0.5 or 1). The explanation lies in the fact 
that cloud Ice Is strongly scattering and in the 
fact that the solar zenith angles are very large 
in the winter polar atmosphere. Without a cloud, 
the solar flux must traverse an effective optical 
depth of r///Q , where r is the vertical optical 
depth measured from the surface to infinity, and 
/iq is the cosine of the solar zenith angle. In 
the winter polar regions, p 0 is quite small (di- 
umally averaged n Q Is approximately 0.32 for 
this case). If a cloud is in the path of the 
solar rays, then most rays are scattered in all 
directions, although with strong forward scatter- 
ing. Some solar rays are scattered to space and 
lost. However, a significant percentage of the 
solar photons are scattered downward at some 
angle 9 which is usually less than the solar 
zenith angle. As a result, r/cos $ is actually 
less than the optical depth traversed had the 
cloud not been there. When solar zenith angles 
are as large as they are in the winter polar 
regions, this is an important effect. This ef- 
fect is even more dramatic farther poleward. The 
same effect Is easily noticed on Earth, especial- 
ly just after sunset. If a cloud- Is present, 
then the surface Illumination can be significant- 
ly greater than If no cloud were present because 
the cloud scatters the solar flux down to the 
surface. 

For low dust opacities, the effect of lowering 
the effective solar zenith angle Is outweighed by 
the scattering of solar flux to space. For mod- 
erate or high dust opacities, lowering the effec- 
tive solar zenith angle decreases the effective 
optical depth of dust that must be traversed and 
allows more solar radiation to reach the sur- 
face. This Is even more Important in NIR and UV 
wavelengths, where C0 2 and 0^ have absorption 
bands. This effect was noted to be particularly 
important for ozone photochemistry on Mars [Lind- 
ner i 1988]. Thermal emission by clouds Increases 
the IR flux absorbed by the polar cap, most no- 
ticeably In the 21- to lOO-pm wavelength range. 
Thermal emission by clouds In the 15-/xm and can 
get absorbed by C0 2 before striking the surface. 

At 70*N latitude, the radiative effect is also 
neutral for small cloud opacities and all dust 
opacities, as shown in Table 5. However, moder- 
ate and large clouds cause up to a 15% increase 
in the flux absorbed by the polar cap. At 70*N 


TABLE 4. Flux Absorbed by the Polar Cap Near Its Edge for Various Cloud and Dust Opacities 

(57*N Latitude. L s - 343*) 


Cloud 


Dust 

Opacity 


0.0 

0.2 

0.5 

1.0 

Opacity 

0.0 

348.177 

340.893 

333.043 

326.734 

0.2 

339.769 

336.545 

333.569 

330.136 

0.5 

333.106 

334.744 

334.143 

332.564 

1.0 

329.832 

333.467 

338.226 

339.172 / 


The vertical optical depth is given. Flux values are given in units of J cm' 2 sol 


Lindner: Martian Polar Cap 


1373 


TABLE 5. Flux Absorbed by the Polar Cap Near the 
Edge of Polar Night for Various Cloud and Dust 
Opacities (70’N Latitude. L s - 343") 


Dust Opacity 


Cloud 


Opacity 

0.0 

0.2 

0.5 

0.0 

175.220 

168.571 

173. A71 

0.2 

167.139 

170.336 

181.032 

0.5 

169. A78 

178.580 

192.287 

1.0 

18A.A56 

208.631 

195.178 


The vertical optical depth is given. Flux 
values are given in units of J cm’ sol’ . 


latitude the solar flux traverses approximately 8 
times as much atmosphere as it vould if it tra- 
veled vertically. The cloud acts to shorten the 
path of the solar radiation, allowing for less 
absorption on the path toward the surface and 
accounting for the increase in absorbed flux. 
This Is especially noticeable in the UV, where 
the flux reaching the surface is 50% higher when 
the radiative effects of a cloud of opacity 1 and 
dust of opacity 1/2 are included. The total 
solar flux which reaches the surface at 70"N lat- 
itude ls only reduced by 10% from "clear” atmo- 
spheric conditions when a cloud opacity of 1 and 
dust opacity of 1/2 are included, which is more 
than offset by the increase in thermal flux, even 
though colder atmospheric temperatures exist at 
70*N latitude. 

It Is at the polar night latitudes that the 
radiative effects of clouds and background dust 
are most significant, as shown in Table 6 (solar 
fluxes are obviously zero). IR fluxes absorbed 
by the polar cap increase up to an order of mag- 
nitude as the cloud opacity and dust opacity 
increase. While dust optical depths of 0.5 are 
shown for completeness, they are unlikely except 
during dust storms. As discussed earlier, dust 
opacities may be lower in the winter polar atmo- 
sphere than over the Viking landers, making back- 
ground dust opacities of 0.5 unlikely. Also, 

TABLE 6 . Flux Absorbed by the Polar Cap in 
Polar Night for Various Cloud and Dust 
Opacities (90*N Latitude. L & - 343*) 


Dust Opacity 


Cloud 


Opacity 

0.0 

0.2 

0.5 

0.0 

14.714 

43.660 

74.765 

0.2 

36.506 

60.984 

88.694 

0.5 

60.951 

81.669 

105.963 

1.0 

90.255 

107.864 

128.389 

3.0 

152.423 

164.712 

178.094 


The vertical optical depth i| given. Flux 
values are given in units of J era sol 


while cloud optical depths as large as 3 could 
occur, optical depths of 0.5 are probably more 
common. However, even simply a cloud optical 
depth of 0.5 and dust optical depth of 0.2 would 
result in the sublimation of the order of A0 cm 
of CO 2 ice over l Martian year. When inserted in 
a polar cap/atraospheric pressure model, this 
would accelerate the predicted loss of all the 
COj Ice on the northern polar cap by approximate- 
ly A0 sols (based on ice depth versus season as 
shown by Leighton and Murray [1966], Cross 
[1971], Briggs [1974], Davies et al. [1977], 
Lindner [1985], and Jakosky and Haberle [this 
issue J) . 

As discussed in the introduction, while polar 
cap/a tmospherlc pressure models have been unable 
to explain the dichotomy of the residual polar 
caps, these models have been fairly successful at 
reproducing observations of atmospheric pressure 
and polar cap recession up to 75" or 80" lati- 
tude. The radiative effects of clouds and back- 
ground dust have the added benefit of not signi- 
ficantly affecting polar cap/atmospheric pressure 
model predictions of the annual variation In 
atmospheric pressure. Clouds and background dust 
have a neutral effect along the cap edge, becom- 
ing more important as one approaches the pole, 
particularly within 80" latitude. The Integral 
of seasonal CO 2 ice over the 80* S to 90* S lati- 
tude band (where the residual polar cap exists) 
Is an order of magnitude less than the integral 
of seasonal C0£ Ice over the planet [Leighton and 
Murray, 1966; Briggs, 1974; Davies et al., 1977; 
Lindner, 1985]. Hence, a process which affects 
primarily the residual polar cap area will not 
significantly affect model predictions of atmo- 
spheric pressure. Also, the rate of recession of 
the seasonal polar cap vould remain essentially 
the same as that predicted by polar cap/a tmo- 
spherlc pressure models equatorvard of 75* lati- 
tude, where good agreement with observations was 
previously obtained. The behavior near the pole 
vould be significantly modified, making for bet- 
ter agreement with observations than previously 
obtained. The rate of recession of the polar cap 
Is an important test, especially since it repeats 
reliably over several years of observation [Ve- 
verka and Coguen, 1973; James et al., 1979; James 
and Lumme, 1982; James, 1982]. 

The Effect of Albedo and Temperature on Qur 

Results 

A polar cap albedo of 0.5 at visible wave- 
lengths was selected for the calculations pre- 
sented here, which ls appropriate during the sub- 
limation phase of the polar cap (late winter, 
spring) [Kleffer, 1979; James and Lumme, 1982]. 
One would expect higher albedos during the forma- 
tion phase of the polar cap (fall, early winter) 
because of the freshly formed ice; unfortunately, 
observations of the polar cap albedo during the 
formation phase of the polar cap are not avail- 
able, due to the extensive cloud cover at this 
time. Higher polar cap albedos would lessen the 
amount of absorbed solar flux, making the thermal 
effects of clouds and dust relatively more Impor- 
tant to the energy budget of the polar cap. A 
polar cap albedo of I would change the pofar cap 
energy budget by nullifying the solar contribu- 
tions in Tables l and 2. For dust-free, cloud- 
free cases, the total flux absorbed by the polar 
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cap at 57’N latitude in late winter would be only 
5% of the value obtained for an albedo of 0.5. 
However , a scenario with a cloud of vertical 
optical depth 1 and dust of vertical optical 
depth 0.5 at 57"N latitude in late winter would 
still result in a flux absorbed by the polar cap 
for an albedo of 1 that is 50% of the flux ab- 
sorbed for an albedo of 0.5, Hence, the presence 
of clouds and dust lessens the ability for a high 
albedo to preserve C0 2 Ice. Furthermore , polar 
cap albedo does not have a linear effect on the 
energy budget, as clouds and dust will reflect 
some of the flux reflected from the polar cap 
back onto the polar cap. Also, a high albedo at 
solar wavelengths is of no consequence at polar 
night latitudes where the Sun never shines. 

Similarly, a low albedo at solar wavelengths 
is not as damaging for the preservation of C0 2 
ice as might be expected, because a substantial 
portion of the energy budget of the polar cap is 
received from IR emission by clouds and back- 
ground dust. Hence, any increase in absorbed 
flux at solar wavelengths from a decrease in 
albedo is not as important to the energy budget 

^ would be in the absence of clouds and 
dust. A polar cap albedo of zero would approxi- 
mately double the absorbed solar flux given in 
the tables, but the percentage Increase to the 
total energy budget is smaller when the effects 
of clouds and dust are included. Furthermore, a 
low albedo at solar wavelengths has no importance 
at polar night latitudes. 

Atmospheric temperatures used in this model 
for the edge of the polar cap were some of the 
coldest observed at polar cap edge latitudes by 
the Viking orbiter [Lindal et al., 1979], Higher 
atmospheric temperatures than those used here 
would actually increase the thermal emission by 
dust and cloud even more, making their radiative 
effects yet more important to the overall energy 
budget. However, the temperatures used are not 
as important as one might expect, because most 
atmospheric dust and cloud lie near the surface 
(about 60% is within 8 km of the surface; one 
atmospheric scale height at 160 K) , and atmo- 
spheric temperatures near the surface remain very 
cl*ose to the C0 2 frost point temperature during 
winter except during dust storms [Lindal et al.. 
1979]. 

Implications for the Polar Cap Asymmetry 

As ozone is more prevalent at northern lati- 
tudes [Barth et al., 1973], ozone was proposed to 
be partly responsible for the hemispherical asym- 
metry In residual polar caps [Kuhn et al. p 
1979]. The maximum ozone abundance observed by 
Mariner 9 is used to test this hypothesis. High- 
er ozone abundances probably exist under cloudy 
and dusty conditions through which the reflec- 
tance spectroscopy technique employed by Mariner 
9 researchers would have had difficulty detecting 
ozone [Lindner, 1988]. However, both ozone heat- 
ing of the atmosphere and the reduction In UV 
light striking the surface would be less Impor- 
tant under very cloudy, dusty conditions [Lind- 
ner, 1988). 

Ozone does absorb most of the ultraviolet 
insolation, which increases the lifetime of C0 ? 
Ice. However, ultraviolet Insolation is only 1% 
of the total Insolation. Thus, even if 0^ ab- 


sorbed all ultraviolet Insolation, C0 2 Ice life- 
times would only be marginally changed. Ozone 
heating of the atmosphere is also minor compared 
to C0 2 and dust heating [Lindner and Thomas, 
1983; Lindner, 1985]. Hence, ozone would not 
cause any appreciable hemispherical asymmetry In 
atmospheric temperature, and therefore In the 
atmospheric thermal emission which strikes the 
polar cap, or in any C0 2 condensation which may 
take place in the atmosphere. 

However, the cloud and dust results have more 
Important Implications regarding the hemispheri- 
cal asymmetry of the residual polar caps. Al- 
though the results are dependent on the polar cap 
albedo and atmospheric temperatures. It does 
appear that the radiative effects of both clouds 
and background dust will appreciably increase 
sublimation of the polar cap at polar night lati- 
tudes. Hence, the radiative effects of clouds 
and background dust would be most noticeable near 
the residual polar cap, which spends all of the 
winter in polar night. As discussed in the in- 
troduction, observational evidence suggests a 
greater abundance of clouds over the northern 
cap than over the southern polar cap. 
Moreover, as also discussed in the introduction. 
It is also expected that the formation and subli- 
mation phases of the seasonal southern polar cap 
will occur with higher background dust opacities 
than the formation and sublimation phases of the 
seasonal northern polar cap. Both of these hemi- 
spherical asymmetries would preferentially in- 
crease the fluxes absorbed by the polar cap at 
polar night latitudes in the north, particularly 
by the residual polar cap area in the north. 

Furthermore, the surface pressure during win- 
ter months is lower over the southern residual 
polar cap than the northern residual polar cap, 
due to differences in elevation [ffoiceshyn, 1974; 
Lindal et al. , 1979] and due to the annual cycle 
in surface pressure [Hess et al., 1980]. If mix- 
ing ratios of cloud and dust are assumed to be 
the same for both poles, then the pressure dif- 
ferential between the poles would result In a 
dust and cloud opacity over the residual polar 
cap In the south that Is 30% less than over the 
residual polar cap in the north during the cru- 
cial winter months. Mixing ratios of dust may 
actually be lower over the southern polar cap 
than over the northern polar cap during winter 
months because the mechanisms for raising and 
maintaining atmospheric dust are pressure depen- 
dent [Pollack et al., 1976, 1979; Pollack and 

Toon, 1982], Lower mixing ratios over the south- 
ern pole would accentuate the asymmetry In cloud 
and dust opacity due to the asymmetry In surface 
pressure. 

The net result of all these hemispherical 
asymmetries is that the solar and infrared flux 
absorbed by the residual polar cap area In the 
north during the winter is greater than that ab- 
sorbed by the residual polar cap area in the 
south. This would allow relatively less frost to 
accumulate over the residual polar cap In the 
north during the fall and winter than in the 
south and cause the C0 2 surface ice to sublimate 
more rapidly in the north during the following 
spring and summer. These suspected hemispherical 
asymmetries in background dust and atmospheric 
pressure have not been included in prior polar 
cap/atraospherlc pressure studies. 
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Other Processes Neglected by Prior Polar 
Cap/Atmospheric Pressure Models 

The radiative effects of cloud and dust are 
not the only processes responsible for the hemi- 
spherical asymmetry in the residual polar caps. 
Prior polar cap/atmospheric pressure studies 
determined that COj ice would completely sublime 
away from both polar caps. This study has shown 
that clouds and background dust serve to Increase 
sublimation over the poles, albeit to a greater 
extent over the northern pole. Some other phe- 
nomenon that has not been included in prior polar 
cap/atmospheric pressure models is causing the 
ice to sublime more slowly, to allow * cc to 
survive southern summer. 

Previous polar cap/atmospheric pressure models 
have assumed that all absorption of sunlight 
occurs at the top surface of the polar cap, when 
in reality it can be distributed over as much as 
several meters depth [Clow, 1987], This simpli- 
fied treatment of ice mlcrophyslcs by prior mo- 
dels yields incorrect sublimation rates for three 
reasons. First, the rate of sublimation/conden- 
sation of a column of Ice depends on the Integral 
of the energy Imbalances at each level within the 
ice, which is not the same as assuming that all 
sources and sinks occur at the top surface. 
Secoud, accurate treatment of Ice microphysics is 
particularly important for thin ice, as much of 
the incoming radiation would go through the cap 
and get absorbed by the surface under the cap . 
Considering that the seasonal cap is at most 1 m 
thick [Leighton and Murray, 1966; Briggs, 1974; 
Davies et al., 1977; Lindner, 1985], this effect 
alone should extend the predicted stability of 
surface ice in the summer. Third, previous mo- 
dels of the sublimation of Martian ice used an 
empirical relationship for sublimation developed 
for turbulent -driven evaporation from terrestrial 
oceans [Clow, 1987]. This parameterization fails 
to account for environmental differences between 
the Earth and Mars and differences between aero- 
dynamlcally rough and smooth surfaces or atmo- 
spheric stability conditions; and the predicted 
sublimation rates are substantially In error (G. 
Clow, personal communication, 1988). 

Furthermore, all prior models Ignored COj 
absorption by the atmosphere in the NIR, which is 
more important at winter polar latitudes where 
the slant paths are so large [Lindner, 1985]. As 
much as 10% of the solar flux can be absorbed, 
which would allow CO 2 Ice to linger longer Into 
the summer months. Moreover, heating within the 
polar cap In near- inf rared wavelengths Is actual- 
ly larger than heating in visible wavelengths in 
the top few millimeters of ice, because near- 
infrared radiation does not penetrate the ice as 
deeply as visible radiation [Clow, 1987]. This 
accentuates the effect that atmospheric absorp- 
tion of near- infrared sunlight has on polar cap 
sublimation. 

Previous models assumed a Lambertian (isotro- 
pic) surface albedo. However, ice does not scat- 
ter isotropically, particularly for solar zenith 
angles above 66* where ice is strongly forward 
scattering [Taylor and Stowe, 1984]. Since the 
solar zenith angle in the winter polar region of 
Mars is usually much greater than 66*, the as- 
sumption of an isotropic albedo for the polar cap 
is in error. This error becomes important when 


the atmospheric feedbacks of clouds and dust 
become important, because clouds and dust do 
scatter a significant amount of solar radiation 
which has been reflected by the polar cap back 
onto the polar cap. This is particularly impor- 
tant for high-albedo objects such as the polar 
caps and the polar hood where solar radiation can 
reflect back and forth many times before escaping 
to space or getting absorbed. 

Furthermore, the total albedo of ice has a 
strong dependence on solar zenith angle [Llan and 
Cess, 1977; Vlscombe and Warren, 1980; Squyres 
and Veverka, 1982; Paige 1985], which previous 
polar cap/atmospheric pressure models have Ignor- 
ed. Ice brightens significantly at high solar 
zenith angles. This Is an Important effect to 
consider in the winter polar regions where solar 
zenith angles vary significantly. The Inclusion 
of the dependence of ice albedo on the solar 
zenith angle should help explain why COj ice is 
retained year-round at the south pole. Since the 
albedo of ice increases and becomes more forward 
scattering at higher solar zenith angles, and 
since the solar zenith angle becomes higher as 
one approaches the pole, the effect will cause 
the albedo to be greatest at the pole. This will 
decrease absorption of sunlight, hence increasing 
survivability of C <>2 ice. Furthermore, at the 
pole the Sun is above the horizon all day, but at 
a very high zenith angle. As one goes toward the 
equator, the Sun spends less time at high zenith 
angles. This will accentuate the solar zenith 
angle effect on ice albedo at the pole even more. 

This is particularly important when viewed 
together with the radiative effects of clouds and 
dust. Assuming that the hemispherical asymme- 
tries in clouds and dust mentioned earlier do 
indeed exist, then the combination of the effects 
of solar zenith angle on albedo and the effects 
of clouds and dust could act to extend the life- 
time of CO 2 ice on the south pole relatively more 
than on the north pole. Another positive aspect 
of solar zenith angle effects on Ice albedo is 
that the extension of CO 2 ice lifetimes at the 
pole due to the inclusion of solar zenith angle 
effects on Ice albedo should not appreciably 
change model predictions of the annual cycle of 
pressure or polar cap recession equatorward of 
75* latitude, since approximately 90% of the sea- 
sonal CO 2 frost Is equatorward of 00* latitude 
[Leighton and Murray, 1966; Briggs, 1974; Davies 
et al., 1977; Lindner, 1985], Hence, the good 
model agreement to polar cap recession equator- 
ward of 80* latitude is retained there. 

Recent studies of the wavelength dependence 
for the albedo of Ice should also be incorporated 
in future polar cap/atmospheric pressure studies, 
given the direct relationship between ice albedo 
and Ice sublimation. Compilations of the optical 
constants of water ice and CO 2 Ice exist from UV 
to microwave wavelengths [Warren, 1984, 1986], 
Not only does experimental work report the wave- 
length dependence of COj/^O Ice mixtures [Klef- 
fer, 1970], but experimental studies also docu- 
ment the change in this wavelength dependence for 
various Ice/dust mixing ratios, various radius 
sizes of Ice particles, and various dust, types 
[Wiscombe and Warren, 1980; Warren and Vlscombe, 
1980; Clark, 1981]. This Is of particular Impor- 
tance In the near- Infrared where the albedo of 
Ice Is highly variable, since significant subll- 
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matlon of the polar cap occurs in the near-infra- 
red [Clow, 1987 ] . The albedo of ice is sensitive 
to grain size in the near- inf rared and sensitive 
to the dust/ice mixing ratio in the visible [War- 
ren and Viscombe, 1980; Uiscombe and Warren, 
1980] . 

Seasonal and latitudinal variations in ice 
albedo occur on the Earth due to changes in dust 
mixing ratios and ice age [Wiscombe and Warren, 
1980; Warren and Wiscombe, 1980 J and should also 
exist on Mars. Latitudinal variations in ice 
albedo have been inferred from Mariner 7 spectra 
[Martin, 1988], and a consistently higher albedo 
for surface ice is inferred to exist in the 
southern hemisphere based on Viking orbiter ob- 
servations [Paige, 1985]. This would lower the 
energy absorbed by the southern polar cap and 
allow the C0£ ice to survive the summer. Season- 
al and latitudinal variations in polar cap albedo 
have been noted to make better agreement between 
polar cap/atmo spheric pressure models and obser- 
vations of atmospheric pressure [Lindner, 1986]. 

Previous polar cap/atmospheric pressure models 
have only examined the condensation of frost and 
have ignored snowfall. Snowfall could be an 
appreciable component of the polar cap [Pollack 
and Haberle, 1988], as it is on Earth, and there- 
fore could account for the lingering of carbon 
dioxide Ice. Carbon dioxide clouds probably 
exist over part of the polar cap [Briggs and 
Leovy, 1974; Hunt et al. , 1980; Paige, 1985], and 
carbon dioxide snow could also be expected to 
form. Furthermore , snowfall should occur mostly 
at polar night latitudes. Sunlit polar winter 
latitudes have a strong inversion which keeps 
atmospheric temperatures well above the C0 2 con- 
densation temperature even while frost is conden- 
sing on the surface [Llndal et al., 1979], How- 
ever, in polar night the heat sources which sup- 
port the temperature inversion at lower latitudes 
are not as effective, allowing for lower atmo- 
spheric temperatures which favor the formation of 
CC >2 clouds and snow* Hence, snowfall might not 
be an important contributor to the polar energy 
budget at latitudes equatorvard of 70*, which is 
where most of the CO 2 frost occurs [Leighton and 
Murray, 1966; Briggs, 1974; Davies et al., 1977; 
Lindner, 1985]. Thus, including snowfall in 
polar cap/atmo spheric pressure models will not 
appreciably change the predictions of the annual 
cycle in pressure or the predictions of polar cap 
recession equatorvard of 70* latitude, both of 
which agree fairly well with observations. How- 
ever, snowfall might be appreciable poleward of 
70* latitude, particularly at 90* latitude, where 
previous models have been incorrect. Snowfall 
would also help counterbalance the radiative 
effects of clouds and dust, which were also most 
important in polar night. Furthermore, southern 
winter is colder and longer than northern winter 
due to the large eccentricity of the Martian 
orbit. Thus, the southern pole might accumulate 
more snowfall than the northern pole, which also 
might help explain why C0 2 frost remains year- 
round at 90*S latitude and not at 90*N latitude. 

Jakosky and Haberle [this Issue] noted that 
prior polar cap/a tmospheric pressure models al- 
ways began with the assumption that the south 
pole was bare of a C0 2 Ice residual polar cap. 
This allowed the surface to warm during summer 
and hence delayed the onset of frost In the 


fall. Jakosky and Haberle showed that the as- 
sumption that a C0 2 ice residual polar cap exist- 
ed in southern summer resulted In C0 2 frost ac- 
cumulating earlier the following fall and hence 
allowed it to survive longer the next summer. 
This idea has the added benefit of being most 
Important at polar night latitudes and would 
therefore help to counter the radiative effects 
of clouds and dust, which are also most important 
at polar night latitudes. 

Recent radar observations of the south polar 
cap suggest the presence of significant roughness 
(D. 0. Muhleman, personal communication through 
B. M. Jakosky, 1988). Theoretically, surface 
roughness In the form of penitentes (spikes of 
ice) may occur in the polar caps [Svitek and Mur- 
ray, 1988]. Surface roughness will change the 
effective solar zenith angle. Surface roughness 
will also allow radiation to be reflected from 
one part of the polar cap directly to another 
part of the polar cap without atmospheric scat- 
tering as an intermediary. . 

Recent calculations of the emissivity of C0 2 
ice as a function of wavelength from the UV to 
the IR under Marslike conditions _shov strong 
variability In emissivity with dust content and 
Ice age (S. Warren, personal communication, 
1988). Under certain conditions, the emissivity 
can be quite low, which would decrease the abili- 
ty of the ice to radiate in the IR. 

Other possibilities include an insulating 
residue on the polar Ice [Saunders et al., 1986], 
which would reduce the rate of sublimation. Wind 
shifting of the polar cap ice may also change the 
recession rate of the polar cap by moving ice 
toward the pole [Briggs, 1974; James et al., 
1979; Kieffer, 1979; Saunders et al., 1985]. 


Summary 

Polar cap/atmospheric pressure models have 
been unable to explain the observed hemispherical 
asymmetry in the composition of the residual 
polar caps on Mars. However, the radiative ef- 
fects of ozone, clouds, and background dust were 
not fully modeled. Ozone has only a minor effect 
on the energy budget of the polar cap. Clouds 
and dust have a net radiative effect that Is 
quite small near the edge of the polar cap, but 
It becomes more substantial in Increasing the 
sublimation rate as one goes poleward. In fact, 
the net radiative . effect of clouds and dust at 
the pole itself increases the loss of C0 2 Ice 
over 1 Martian year of the order of 40 cm or 
more, which would result In the loss of C0 2 ice 
of the order of 40 sols earlier than previously 
predicted . 

However, several other processes have also not 
been fully included in polar cap/atmospheric 
pressure models. These processes Include pene- 
tration of radiation into and through the Ice, 
bidirectional surface reflectance. Ice albedo 
dependence on solar zenith angle, season, and 
latitude, lower Ice emissivity, decreased heat 
conduction from the surface due to existing ice, 
snowfall, surface residues. surface roughness, 
and wind shifting of ice. The importance of each 
of these processes to the energy budget Is cur- 
rently being assessed by B. L. Lindner and B. M. 
Jakosky- Most of these processes would lengthen 
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C0 2 Ice survivablLity either by increasing con- 
densation or decreasing sublimation of C 0 2 ice 
and would counter the increased sublimation caus- 
ed by clouds and dust. However, most of these 
processes work equally effectively at both poles 
and do not generate an asymmetry. Hemispherical 
asymmetries in cloud and dust opacity and in sur- 
face elevation would result in a hemispherical 
asymmetry in the net radiative flux from clouds 
and dust which is absorbed by the polar cap. The 
key point is that the net radiative effect of 
clouds and background dust is to shorten ice 
lifetimes on the north pole relative to the south 
pole, helping to explain the asymmetry in the 
residual polar caps on Mars. Furthermore, the 
radiative effects of clouds and dust do not ap- 
pear to significantly affect model predictions of 
polar cap recession equatorvard of 75 * latitude 
or model predictions of the annual cycle in atmo- 
spheric pressure, both of which agreed fairly 
well with observations. The presence of clouds 
and dust also lessens the ability for a higher 
albedo to preserve CO 5 ice and for a lower albedo 
to increase sublimation, decreasing the impor- 
tance of albedo in determining ice sublimation 
rates . 
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Penetration of Light I nt o the Martian Polar _Cap: 
Implications for the Energy Budget 

B.L. Lindner (AER) and B.M. Jakosky (U. Colo.) 

We have modified the Jakosky and Haberle ( J. Geophys . 
Res . . 95, 1359, 1990) model of the seasonal polar cap 
on Mars to include penetration of solar radiation into 
the cap itself, based on the theoretical work of Clow 
(Icarus, 72, 95, 1987). We find that the inclusion of 
light penetration slightly decreases the albedo needed 
in the model to keep (X^-ice year-round at the south 
pole by on the order of 1%. The required albedo is 
decreased because some solar radiation Is used to heat 
the subsurface, and not all of this heat Is transport- 
ed back to the surface. Furthermore, the number of 
sols for which the north pole Is free of CC^-Ice is 
increased slightly by approx. 5 sols. By allowing 
solar radiation to be absorbed at depth, the surface 
does not become as warm. This decreases the Infrared 
radiation emitted by the surface and increases the 
heat conduction from the subsurface, which in turn 
delays the onset of frost formation in the fall. 

Overall, we conclude that the penetration of light 
into the polar cap has only a small effect on the 
energy budget of the polar cap. Given the 
uncertainties currently present In albedo and other 
parameters, the effect of light penetration is second 
order, and can be neglected In models of the polar cap 
energy budget. 


AN EFFICIENT AND ACCURATE TECHNIQUE TO COMPUTE THE ABSORPTION. EMISSION. AND TRANSMISSION 
OF RADIATION BY THE MARTIAN ATMOSPHERE. Bernhard Lee Lindner. AER Inc.. 840 Memorial Orrve, 
cLSJe MA 02139; Tho^s P. Ackerman. Dope. of Meteor., ^nn. State Unrv. . Unrversrty 
Park. PA 16802; and James B. Pollack, NASA/ARC. Moffett Field, CA 9403b. 

INTRODUCTION. COo comprises 95% of the composition of the martian atmosphere jLJ. However 
the Xia^ ^osphlre .'iso has a high aerosol content. Dust opacities vary fro. less than 
0.2 to greater than 3.0. primarily on a seasonal basis with the occurrence of global dust 
storms during southern spring (2]. Ice-cloud opacities vary from 0 to •*“ l ' ™ 

large amounts occurring at winter polar latitudes (3J . <*>2 Is an act ve a sor er an 

in near-IR and IR wavelengths; the near-IR absorption bands of CO, provide significant h g 

of the atmosphere, and the 15 pm band provides rapid cooling(4-5 J . However dust and ice- 
cloud aerosols have high scattering albedoes in solar wavelengths and are highly absorbing a 
infrared wavelengths, and are as important as C0 2 in the atmospheric energy budget l 5 J- 

Including both CO, and aerosol radiative transfer simultaneously in a model is diffic . 
Aerosol radiative transfer requires a multiple-scattering code, while C0 2 radiative transfer 
must deal with complex wavelength structure, as shown in Fig. 1. The problem can be so ve 
exactly by inserting the CO, absorptance for each spectral line into a multiple-scattering 
code, but the 15,® band alo 2 ne has on the order of 10.000 lines, making such * '^Stole 
tedious and expensive. It is this difficulty of simultaneously treating aerosol multiple 
scattering and the banded absorption structure of 
CO* that prompts most radiative- transfer studies 
of the martian atmosphere to consider either a 
pure-C0 2 or pure -dust atmosphere. this approxi- 
mation simplifies treatment, but is inaccurate. 

One alternative technique that has recently 
been developed for atmospheric applications is 
the exponential -sum or k-distribution approxima- 
tion (8-21] . The transmission of a homogeneous 
atmosphere is actually independent of the order- 
ing of the absorption coefficient, k, in frequen- 
cy space within a spectral interval, depending 
only upon the percentage of the spectral interval 
that has a particular value of k. The percentage 
of the spectral interval which has values between 
k and k + Ak can be formulated in a probability 
density function f(k) shown schematically in Fig. 

2. The chief advantage of the exponential -sum 
approach is that the integration over k space of 
f(k) can be computed more quickly than the inte- 
gration of k^ over frequency. The exponential - 
sum approach is superior to the photon- 
path -distribution and emissivity tech- 
niques for dusty conditions (22, 19, 

23]. Our work is the first application 
of the exponential -sum approach to mar- 
tian conditions. 

THEORETICAL APPROACH . The transmit- 
tance of th e 15 fun band and the near - IR 
bands of C0 2 was computed using the 
FASCOD line-by-line transmittance model 
[24], modified for martian condi- 
tions. Computations with the modified 
FASCOD model were made at 3 tempera- 
tures (125K, 200K, 300K) and 5 pres- 
sures (100 mb, 10 mb, 1 mb, 0.1 mb, 

0.01 mb); these cover the range of tem- 
perature and pressure currently observ- 
ed in the atmosphere at all latitudes, 
seasons, and altitudes up to 40 kra, and 
also can be used for early Mars, dense 
atmosphere studies. The near-IR and 
15pm bands were broken Into spectral 
sub- intervals (see Fig. 1), and the 
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Figure 1. C0 2 15 fim band transmission 
at 20 km altitude looking upward, at 
temperatures of 200 and 300K [25]. 
Also shown are the sub- intervals used 
for the 8 term and 16 term fits. 



Figure 2. A schematic illustration portraying 
the essence of the exponential - sum approach. 

(a) shows a schematic of absorption line spectra 
at two different pressures. In (b) the two pro- 
bability density functions f(k) associated with 
(a) are illustrated. The shaded area depicts 
the strongest absorption (i.e., largest k) for 
the same spectral interval (i.e., f(k) for dif- 
ferent pressures are correlated) . Integration 
of f(k) over k replaces the Integration of k^ 
over v (modified from [28], [29]). 
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transmittance lor each layer, Tr, as a function of C0 2 column abundance within that layer, u, 
was, fit with a series of weighting coef f ic ients , a^, and exponential coefficients, b^z 
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We tried fitting procedures based on Wlscotnbe and Evans (14] and an Improved version of Acker- 
man et al. (13], and found both yielded similar results. Both of these procedures avoid the 
ill-condltlonlng of earlier exponent ial-sum routines, and produce more accurate and unique so- 
lutions (13,14], The exponential -sum fit reproduced the FASOOD transmittances to better than 
10~* for all C0 2 abundances considered. The frequency sub - intervals are picked to try to 
minimize the variation in line strengths within the sub -interval. As shown In Fig. 1, 2 sub- 
intervals covered the band center, 2 covered the far line wings, 2 covered the near line 
wings, and 2 covered the transition from wings to band center. The vertical inhomogeneity of 
the atmosphere Is treated by using homogeneous layers, with the absorption coefficients k for 
all layers correlated in frequency space, i.e. the sub- intervals of the spectral band which 
have the maximum absorption also have the largest k values (Ackerman et al. , 1976). An inter- 
polation is used for temperatures and pressures which fall in between the values at which the 
exponential-sum coefficients were computed. We compared a logarithmic interpolation to a 
linear Interpolation, and found the logarithmic interpolation to be more accurate. 

INCORPORATION IN MULTIPLE -SCATTERING MODELS . Vertical optical depths of C0 2 absorption for 
each term number, frequency sub-interval, and atmospheric layer are bjU. C0 2 is combined with 
dust and cloud in that the total optical depth Tj, single -scattering albedo «£, and phase 
function P for eac^ terjj, frgquengy su^-interval, and layer are given by{16]: 


T. - b.u + r + t 4* t + t + t 
1 1 s a s a 

“1 - <' R + \ + 's> / T i 
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+ r C P C + r¥ 
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where: r R - Rayleigh scattering optical depth for that layer 

- Dust (D) and Cloud (C) scattering (s) and absorption (a) optical depth for 
* that layer 

p D,C,R M dust doud (C), and Rayleigh scattering (R) phase function 

The multiple -scattering code Is run once for each terra In the sura using the T, and P 

appropriate to that term, and the resultant fluxes, F* , (or Intensities) are then summed and 
weighted by a^ to give the total flux, F, (or Intensity; over the frequency sub-interval: 


F - Z a F (b u) 
i-1 


(5) 


Tables of the exponential -sum coefficients, a^ and bj_, 
can be obtained from the authors. 

NUMERICAL STUDIES . The number of terras in the 
series of exponentials (n in equation 1) and the num- 
ber of sub- intervals into which the spectral band is 
broken can be varied to increase the desired accur- 
acy. Figure 3 shows a comparison between the 15 /an 
cooling rates computed using various numbers of fre- 
quency sub- intervals and terms within each sub- inter- 
val. The temperature of the atmospheric model rose 
linearly from 150K at the surface to 160K at 10km 
altitude, and then fell linearly to 130K at 40km alti- 
tude. These temperatures are typical of the winter 
polar atmosphere on Mars. Using 4 terms and 4 sub- 
intervals results in an error of less than 10% In the 
lowest 10 km, but results in substantial errors at 
higher altitudes. The explanation is that the 4 term 
fit could not capture the effects of both the band 
center and the wings at all altitudes, and in this 
case we emphasized the wings, which are more important 
at lower altitudes where the band center is saturated 
(see Fig. I). Adding additional terras from 8 to 16 Is 
far less noticeable. 



Figure 3. 15 ptm band cooling rates 

calculated with an exponential - sum 
approach using 4 frequency sub- In- 
tervals and 4 terms in the sum, 10 
sub - intervals and 8 terms, and 10 
sub- intervals and 16 terms. A win- 
ter polar temperature profile is 
used (see text) , but only 200K expo- 
nential-sura coefficients are used 
for comparison purposes. 
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We also separately derived exponent Lai - sura coeffi- 
cients based on the Cal 'tsev and Osipov [25] line -by- 
line calculations of the 15 fita band. The number of 
terms In the fit (n in equation 1) were varied, to 
confirm our FASCOD work on the number of terras requir- 
ed to give an accurate fit. Figure A shows cooling 
rates computed with several fits to the Gal' tsev and 
Osipov parameter izat ions , as well as to the FASCOD 
transmlttances . Cooling rates compare favorably in 
the lower atmosphere (below 10 km altitude) . Cal tsev 
and Osipov only considered temperatures as cold as 
200K. We extrapolated the temperature dependence to 
colder temperatures [26], Discrepancies between the 
Gal * tsev and Osipov and FASCOD cooling rates are due 
to the inaccuracy in the temperature extrapolation, 
particularly at the higher, colder altitudes. Again, 
note that the A and 5 term fits become inaccurate 
above 10 km altitude, as in Fig. 3. Also, the 15 term 
fit yields no marked Improvement over the 8 term fit, 
as also shown with the FASCOD fits in Fig. 3. 

Cooling rates computed using the FASCOD exponen- 
tial-sum transmlttances also compared well below 10 km 
altitude with cooling rates computed using the Pollack 
et al. [27, A] parame ter Izat ions of CO 2 transmit- 



Coo ling Rate CK/day) 

Figure A. 15 /xm band cooling rates 
calculated with A types of exponen- 
tial -sum coefficients: A and 15 

term fits to the transmission para- 
meter izat ions of [25] (abbrev. C&O) , 
and A and 8 term fits to the FASCOD 
transmission model. The winter po- 
lar atmosphere model is used, with 
logarithmic temperature interpola- 
tion for the exponential coeffi- 
cients . 


tance. At higher altitudes, discrepancies exist due 
to the use of the strong- line approximation by Pollack 
et al, which emphasizes the effects of the line 
wings. However, comparisons of our exponential -sum 

transmlttances with techniques commonly used for the , . 

terrestrial atmosphere have indicated that our exponential -sum transmlttances cou1 ^ b 
error above 20 to altitude. J. Pollack Is seeking to modify our exponential -sum approach to 
use two sets of sums, one applied to the line centers and one applied to the line wings. 
Initial tests seem to show better agreement In the upper atmosphere of Mars. 

B.L. Lindner acknowledges support by NASA contract NASW-4444. 
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SOLAR AND IR RADIATION NEAR THE MARTIAN SURFACE: A PARAMETERIZATION FOR C0 2 TRANSMITTANCE. 
BeiTtihard Lee Lindner, AER Inc., 640 Memorial Drive, Cambridge, MA 02139; Thomas P. Ackerman, Dept, 
of Meteor.. Penn. State Unlv. . University Park, PA 16802; James B, Pollack and 0. Brian Toon, NASA- 
/ARC, Moffett Field. CA 94035; Gary E. Thomas, APAS Dept., Unlv. of Colorado, Boulder, CO 80309 


<X >2 comprises 951 of the atmospheric composition of the martian atmosphere (lj. 
CO 2 Is an active absorber and emitter In near-IR and IR wavelengths; the near-IR absorption bands of 
CO 2 provide significant heating of the atmosphere, and the 15 /jib band provides rapid cooling (2- 
5], However, the Mars atmosphere also has a high aerosol content. Dust opacities vary from less 
than 0.2 to greater than 3.0, primarily on a seasonal basis with the occurrence of global dust 
storms during southern spring (6], Ice-cloud opacities vary from 0 Co greater than 1, with large 
amounts occurring at winter polar latitudes (7J. Dust and Ice-cloud aerosols have high scattering 
albedoes In solar wavelengths, and are highly absorbing at Infrared wavelengths, and are as 
Important as COj In the atmospheric energy budget (5). 

Including both C0j and aerosol radiative transfer simultaneously in a model Is difficult. 
Aerosol radiative transfer requires a multiple -scattering code, while CO 2 radiative transfer must 
deal with complex wavelength structure. The problem can be solved exactly by Inserting the CO 2 
absorptance for each spectral line into a multiple -scattering code, but the 15/ua band alone has on 
the order of 10,000 lines, making such a computation tedious and expensive. It Is this difficulty 
of simultaneously treating aerosol multiple scattering and the banded absorption structure of CO 2 
that prompts most radiative- transfer studies of the martian atmosphere to consider either a pure -CO 2 
or pure-dust atmosphere. This approximation simplifies treatment, but is inaccurate. 

One alternative technique that has recently been developed for atmospheric applications is the 
exponential- sum or k- distribution approximation [8-21]. The transmission of a homogeneous atmo- 
sphere is actually independent of the ordering of the absorption coefficient, k, in frequency space 
within a spectral interval, depending only upon the percentage of the spectral interval that has a 
particular value of k. The percentage of the spectral interval which has values between k and k + 
AT can be formulated in a probability density function f(k) shown schematically in Figure 1. The 
chief advantage of the exponential- sum approach is that the integration over k space of f(k) can be 
c<nnputed more quickly than the integration of k v over frequency. The exponential -sum approach is 
superior to the photon-path-distribution and emissivity techniques for dusty conditions [22, 19, 
23]. Our work is the first application of the exponential -sum approach to martian conditions. 

THEORETICAL APPROACH . The transmit- 
tance of the 15 fm band and the near-IR 
bands of CO 2 was computed using the 
FASCOD line-by-line transmittance model 
[24], modified for martian conditions. 

Computations with the modified FASCOD 
model were made at 3 temperatures (125K, 

200K, 300K) and 5 pressures (100 mb, 

10 mb, 1 mb, 0.1 mb, 0.01 mb); these 
cover the range of temperature and 
pressure currently observed in the atmo- 
sphere at all latitudes, seasons, and 
altitudes up to 40 km, and also can be 
used for early Kars, dense atmosphere 
studies. The near-IR and 15 /jm bands 
were broken into spectral sub- intervals, 
and the transmittance for each layer, 

Tr, as a function of CO 2 column abun- 
dance within that layer, u, was fit with 
a series of weighting coefficients, a^. 



and exponential coefficients, b ^ : Tr(u) 


Figure 1. A schematic Illustration portraying the 
essence of the exponential -sum approach. (a) shows 
a schematic of absorption line spectra at two dif-' 
ferent pressures. In (b) the two probability densi- 
ty functions f(k) associated with (a) are illus- 
trated. The shaded area depicts the strongest ab- 
sorption (l.e., largest k) for the same spectral 
interval (i.e., f(k) for different pressures are 
correlated). Integration of f(k) over k replaces 

the integration of k over u (reproduced from [28]). 


sum(a i exp(-b^u) ) from term number 
1-1, n. We tried fitting procedures based on Uiscombe and Evans [14] and an Improved version of 
Ackerman et al. [13], and found both yielded similar results. Both of these procedures avoid the 
Hl-condltloning of earlier exponential -sum routines, and produce more accurate and unique ^solutions 
[13,14]. The exponential -sum fit reproduced the FASCOD transmittances to better than 10 for all 
COj abundances considered. The sub- intervals are picked to try to minimize the variation in line 
strengths within the sub- interval. Vertical optical depths of CO 2 absorption for each term number, 
frequency sub- Interval , and atmospheric layer are b^u. This CO 2 absorption opacity and the Rayleigh 
scattering opacity Is combined with dust and cloud opacities, single- scattering albedoes and phase 
functions (e.g. t [16]). The multiple -scattering code is run once for each term In the sum using the 
T, 0 , and P appropriate to that term, and the resultant fluxes (or Intensities) are then summed and 
weighted by the to give the total flux (or intensity) over the spectral sub - Interval . An 

Interpolation Is used for temperatures and pressures which fall in between the values at which the 
exponential -sum coefficients were computed. We compared a logarithmic Interpolation to a linear 
interpolation. and found the logarithmic Interpolation to be more accurate. Tables of the 
exponent Ial - sum coefficients, a^ and b^, can be obtained from the authors. 
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C0 2 TRANSMITTANCE IN THE MARTIAN ATMOSPHERE: 


RESULTS . The number of terns in the series of expo- 
nentials and the number of sub- intervals into which the 
sp'ectral band is broken can be varied to increase the 
desired accuracy. Figure 2 shows a comparison between 
the 15 jam cooling rates computed using various numbers 
of frequency sub- intervals and terms within each sub- 
interval. The temperature of the atmospheric model 
rose linearly from 150K at the surface to 160K at 10km 
altitude . and then fell linearly to 130K at 40km alti- 
tude. These temperatures are typical of the winter 
polar atmosphere on Mars. Using 4 terms and 4 sub- 
intervals results in an error of less than 10Z in the 
lowest 10 km, but results in substantial errors at 
higher altitudes. The explanation is that the 4 term 
fit could not capture the effects of both the band cen- 
ter and the wings at all altitudes, and in this case we 
emphasized the wings, which are more important at lower 
altitudes where the band center Is saturated. Adding 
additional terms from 8 to 16 is far less noticeable. 
However, very recent work by J. Pollack comparing these 
exponential -sum transmit tances with techniques commonly 
used for the terrestrial atmosphere have indicated that 
our exponential -sum transmittances could be In error 
above 20 km altitude. J. Pollack Is seeking to modify 
our exponential- sum approach to use two sets of sums , 
one applied to the line centers and one applied to the 
line wings. Initial tests seem to show better agree- 
ment in the upper atmosphere of Mars. 

We also separately derived exponential -sum coeffi- 
cients based on the Gal'tsev and Osipov (25] Itne-by- 
llne calculations of the 15 jim band. The number of 
terms In the fit were varied, to confirm our FAS COD 
work on the number of terms required to give an accu- 
rate fit. Figure 3 shows cooling rates computed with 
several fits to the Cal'tsev and Osipov parameteriza- 
tlons, as veil as to the FAS COD transmittances. Cool- 
ing rates compare favorably In the lower atmosphere 
(below 10 km altitude). Gal'tsev and Osipov only con- 
sidered temperatures as cold as 200K. We extrapolated 
the temperature dependence to colder temperatures 
[26]. Discrepancies between the Gal'tsev and Osipov 
and FASCOD cooling rates are due to the Inaccuracy in 
the temperature extrapolation, particularly at the 
higher, colder altitudes. Again, note that the 4 and 5 
term fits become Inaccurate above 10 km altitude, as In 
Figure 2. Also, the 15 term fit yields no marked 
Improvement over the 8 term fit, as also shown with the 
FASCOD fits in Figure 2. 

Cooling rates computed using the FASCOD exponential - 
sum transmittances also compared well below 10 km alti- 
tude with cooling rates computed using the Pollack et 
al. [27,4] parameterizations of CO 2 transmittance. At 
higher altitudes, discrepancies exist due to the use of 
the strong- line approximation by Pollack et al, which 
emphasizes the effects of the line wings. 



Cooling Rate (*K/day) 

Figure 2. 15 band cooling rates cal- 

culated with an exponential -sum approach 
using 4 frequency sub -Intervals and 4 
terms in the sum, 10 sub- intervals and 8 
terms, and 10 sub- Intervals and 16 
terms. A winter polar temperature pro- 
file Is used (see text) , but only 200K. 
exponential -sum coefficients are used 
for comparison purposes. 



Cooling Rate ( *K/day) 
Figure 3. 15 /im band cooling races 

calculated with 4 types of exponen- 
tial-sum coefficients: 4 and 15 term 

fits to the transmission parameteriza- 
tlons of [25] (abbrev. G&O) , and 4 and 
8 term fits to the FASCOD transmission 
model. The winter polar atmosphere 
model Is used, with logarithmic tern-, 
perature interpolation for the expo- 
nential coefficients. 

B.L. Lindner acknowledges support 
by NASA contract NASV-4444. 
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C0 2 TRANSMITTANCE IN THE MARTIAN ATMOSPHERE: 
AN EXPONENTIAL- SUM FIT FOR USE IN MULTIPLE- 
SCATTERING MODELS 

Bernhard Lee Lindner (Atmospheric and Environ- 
mental Research Inc. , 840 Memorial Drive, 

Cambridge, MA 02139, U.S.A.) 

James B. Pollack (NASA/Ames Research Center, 
Moffett Field, CA 94035, U.S.A.) 

Thomas P. Ackerman (Dept, of Meteorology, Penn. 
State Univ. , University Park, PA 16802 USA) 

C0 2 and airborne dust are both important in 
heating and cooling the martian atmosphere. 
However, theoretical modeling of both C0 2 and 
dust radiative transfer simultaneously is dif- 
ficult because dust radiative transfer requires 
a multiple scattering code, while C0 2 radiative 
transfer must deal with complex wavelength 
structure. 

We have approximated the C0 2 transmittance at 
near-IR and IR wavelengths In the lower martian 
atmosphere with an exponential -sum fit. Expo- 
nential-sum fitting allows for the incorporation 
of C0 2 absorption and emission In a multiple 
scattering computer model in a straightforward, 
efficient and accurate manner. This makes it 
possible for C0 2 and dust heating and cooling of 
the martian atmosphere to be easily treated 
s imul taneously . 

Comparison of our C0 2 cooling and heating 
rates to those derived from several other tech- 
niques shows good agreement in the lower atmo- 
sphere below 20 km altitude , particularly near 
the surface. An improved version of the expo- 
nential-sum approach is being developed to treat 
higher altitudes as well. 

Tables of exponential -sum coefficients are 
vailable from the authors. 
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MARS SEASONAL C0 2 -ICE LIFETIMES AND THE ANGULAR DEPENDENCE OF ALBEDO 
Bernhard Lee Lindner . AER, 840 Memorial Drive, Cambridge MA 02139 USA 

The albedo of the polar caps on Mars brightens appreciably at high 
solar zenith angle (Warren et al. f J . Geophys . Res , . 95, 14717, 1990), 
an effect not included in prior polar-cap energy-balance models. This 
decreases absorption of sunlight by the polar cap, hence decreasing 
sublimation of C0 2 ice. Lindner ( J . Geophys . Res . « 95 . 1367, 1990) 
has shown that the radiative effects of clouds and airborne dust will 
increase sublimation of C0 2 ice over that predicted by prior polar-cap 
energy-balance models. Furthermore, observations hint that more 
clouds may exist in the northern hemisphere, which Lindner (1990) has 
shown would sublime C0 2 ice more quickly in the north than in the 
south. I show here that the effects of the solar zenith angle depen- 
dence of albedo and the radiative effects of clouds and dust offset 
each other, but act to extend the lifetime of C0 2 ice on the south 
pole more than on the north pole, possibly explaining the observed 
hemispherical asymmetry in the residual polar caps without the need of 
a hemispherical asymmetry in polar-cap albedo required by prior 
models. Another positive aspect of this solution is that neither the 
Inclusion of the solar zenith angle dependence of albedo nor the radi- 
ative effects of clouds and dust should appreciably change prior model 
agreement with observations of the annual cycle of surface pressure 
and the recession of the polar caps equatorward of 75° latitude. 
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WHY IS THE NORTH POLAR CAP ON MARS DIFFERENT THAN THE SOUTH POLAR CAP? 

Bernhard Lee Lindner . Atmospheric and Environmental Research, Inc. 

840 Memorial Drive, Cambridge, Mass. 02139-3794, USA 

Introduct Ion . One of the most puzzling mysteries about the planet Mars is 
the hemispherical asymmetry In the polar caps. Every spring the seasonal polar 
cap of CO 2 recedes until the end of summer, when only a small part, the residu- 
al polar cap, remains. During the year that Viking observed Mars, the residual 
polar cap was composed of water ice in the northern hemisphere [Kieffer et al., 
Science . 194 . 1341, 1976] but was primarily carbon dioxide ice In the southern 
hemisphere [Kieffer, J . Geoohvs . Res . . £4, 8263, 1979]. Scientists have sought 
to explain this asymmetry by modeling observations of the latitudinal recession 
of the polar cap and seasonal variations in atmospheric pressure (since the 
seasonal polar caps are primarily frozen atmosphere, they are directly related 
to changes In atmospheric mass) . These models reproduce most aspects of the 
observed annual variation In atmospheric pressure fairly accurately. Further- 
more, the predicted latitudinal recession of the northern polar cap in the 
spring agrees well with observations, including the fact that the CO 2 Ice is 
predicted to completely sublime away. However, these models all predict that 
the carbon dioxide ice will also sublime away during the summer In the southern 
hemisphere, unlike what is observed. This paper will show how the radiative 
effects of ozone, clouds, and airborne dust, light penetration into and through 
the polar cap, and the dependence of albedo on solar zenith angle affect CO 2 
Ice formation and sublimation, and how they help explain the hemispherical 
asymmetry In the residual polar caps. These effects have not been studied with 
prior polar cap models. 

Ozone. Clouds, and Airborne Dust . Since 0^ is more prevalent in the 
northern hemisphere than In the southern hemisphere, 0^ was suggested as a 
cause for the hemispherical asymmetry In the residual polar caps by Kuhn et al. 
( J . Geonhvs . Res . . 84 . 8341, 1979). However, Lindner (submitted to Icarus , 
1991) has shown that O 3 has a minor effect on the atmospheric temperature, and 
hence on the infrared radiation which strikes the polar cap, and Lindner (J^ 
Geophvs . Res. . 95, 1367, 1990) has shown that 0 3 absorbs less than 1% of the 
total solar radiation absorbed by the polar cap. Thus, 0 3 Is not an Important 
consideration in the polar cap energy budget. 

Lindner (1990) has computed the solar and thermal flux striking the polar 
cap of Mars for various ozone, dust, and cloud abundances and for three solar 
zenith angles. These calculations have been inserted in the polar-cap models 



of Lindner ( Eos Trans. AGU . 67, 1078, 1986) and Jakosky and Haberle ( J, Geo- 
phvs . Res. . 95, 1359, 1990). Vertical optical depths of dust and cloud ranging 
from zero to 1 cause little change in the total flux absorbed by the polar cap 
near its edge but increase the absorbed flux significantly as one travels pole- 
ward. Observed hemispherical asymmetries in dust abundance, cloud cover, and 
surface pressure combine to cause a significant hemispherical asymmetry in the 
total flux absorbed by the residual polar caps, which helps to explain the 
dichotomy in the residual polar caps on Mars. 

Light Penetration . Penetration of solar radiation into the cap itself is 
included in my polar cap model, based on the theoretical work of Clow (Icarus . 
72 . 95, 1987). I find that the inclusion of light penetration slightly 
decreases the albedo needed in the model to keep CC^-ice year-round at the 
south pole by on the order of 1%. The required albedo is decreased because 
some solar radiation is used to heat the subsurface, and not all of this heat 
Is transported back to the surface. Overall, I conclude that penetration of 
light into the polar cap has only a small effect r on the polarcap energy budget. 

Albedo and the Solar Zenith Angle . Warren et al. ( J . Geonhvs . Res . . 95, 
14717, 1990) has computed the dependence of the albedo of the martian polar 
caps on solar zenith angle, and these calculations have been included In my po- 
lar cap model. Since the albedo of Ice increases and becomes more forward 
scattering at higher solar zenith angles, and since the solar zenith angle be- 
comes higher as one approaches the pole, the albedo is greatest at the pole. 
This decreases absorption of sunlight, hence increasing survivability of CC >2 
ice. In fact, this increases the survivability of ice enough to offset the de- 
crease in survivability of Ice due to the radiative effects of clouds and dust. 

Discussion. The combination of the effects of solar zenith angle on albe- 
do and the radiative effects of clouds and dust act to extend the lifetime of 
CO 2 ice on the south pole relatively more than on the north pole, explaining 
the hemispherical asymmetry In the residual polar caps without the need of a 
hemispherical asymmetry In polar cap albedo. Another positive aspect this 
solution is that neither the inclusion of solar zenith angle effects on Ice 
albedo nor the radiative effects of clouds and dust should appreciably change 
model predictions of the annual cycle of pressure or polar cap recession equa- 
torward of 75° latitude, since approximately 90% of the seasonal CO 2 frost is 
equatorward of 80° latitude. Hence, the good model agreement noted by prior 
researchers to the seasonal cycle in atmospheric pressure and to the recession 
of the polar cap equatorward of 80° latitude is retained. 



